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^ ■ Abstract 
• I— I . 

^ ' There are two schools, or hnes of thought, that try to unify the 

■ apparently divergent laws of dynamics and thermodynamics and to 

explain the observed time-asymmetry of the universe, and most of its 
sub-systems, in spite of the fact that these systems are driven by time- 
symmetric evolution equations. They will be called the coarse-graining 
and the fine-graining schools (even if these names describe only a part 
of their philosophy). Coarse-graining school obtains time-asymmetry 
via a projection of the state space on a space of "relevant" states. 
The corresponding projection of the primitive reversible evolution laws 
yields effective irreversible evolution laws for the relevant states. Fine- 
graining always use the same primitive reversible evolution laws. But 
these laws (in adequate extensions of the usual spaces where these laws 
are formulated) have a set of solutions S that can be decompose in 
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two subsets 5+ and S- of time asymmetric solutions. Choosing one of 
these two sets, as the arena to formulate the theory, time asymmetry 
is established. The aim of these lectures is to explain, in the simplest- 
self-contained, unbiased, and, honest way, the main characteristics of 
both schools and to point out the advantages and disadvantages of 
both formalism, in such a way that, the polemic between the schools, 
turns out to be explicit and organized in the mind of the reader (who 
will be considered the supreme judge to give the final verdict). 

Some cosmological features of the theory will be also considered, 
mainly the problem of the low entropy initial state of the universe 

Dynamics, thermodynamics, and time-asymmetry 



• PACS Nrs. 05.20-y, 03.65, BZ, 05.30-d 



1 Introduction. 

In these lectures we will study and try to solve two, long standing, problems 
of theoretical physics. 

1.1 The problem of time asymmetry. 

The problem of the existence of the arrow of time or, what is the same 
thing, the problem of time asymmetry of the universe, can be stated in two 
questions: 
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i. How can it be that the universe is time- asymmetric if all the relevant 
physical laws are time symmetric? 

ii. Why all arrows of time point in the same direction? 

In fact, universe has several time asymmetries, namely the various arrows 
of time: thermodynamical, electromagnetical, psicological, etc., while the 
main laws of nature are timc-symmctric (because, as usual, wc will neglect 
the laws of weak interaction, since it is very difficult to imagine a mechanism 
that explain the time asymmetry of the universe based in these laws [1]). 

In these lectures we would like to answer these questions given an ad- 
equate mathematical formalism to the problem and using several, old and 
new, well known ideas ([2], [3], [4]). In doing so we must ffist precise two 
important words: conventional and substantial ([1],[5]). In mathematics we 
use to work with identical objects, like points, the two directions of an axis, 
the two semi cones of a null cone, etc. In physics there are also identical 
objects: like identical particles, spin directions, etc. -Among identical objects 
there is always a mathematical transformation that exchange these object 
leaving the system unmodified. If we are forced to call identical objects with 
different names we will say that be are establishing a conventional difference 
among these objects, e. g. when we call -|- and - the two directions of an axis 
or" past" and "future" the two semi cones of a null cone. If some objects are 
different we will say that there is a substantial difference among them. The 
problem of the arrow of time is that past and future are only conventionally 
different, in usual physical theories, while we have the filling that past and 
future are substantially difference, in facts in the past events had happened, 
while in the future event could only happen. 

In theories endowed with time-symmetric evolution equations, as those 
we will deal with, it is quite impossible to find time substantial asymmetry 
using rigorous mathematical manipulations. But normally we can find in 
these theories, as we shall see, two identical mathematical structures, one 
related with the past and one related with the future, e. g. two subspaces of 
the space of solutions of the theory. Nevertheless these structures are only 
conventionally different, because they are related with a time inversion. But 
within these structures past is substantially different than future. To chose 
one of these structures,or the other, is physically irrelevant, since time inver- 
sion exchange one structure with the other, leaving the universe unchanged. 
Therefore to create an arrow of time we just conventionally chose one of the 
structures. This choice is irrelevant, as irrelevant as to chose one face of a 
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dice if all faces are marked with the same number. But when we have chose 
one of the structures a substantial difference is also created, between past 
and future, within these structure,and an arrows of time appears. This is the 
method we will use to create all the arrows of time, both in coarse-graining 
and fine-graining cases (see section 6). 

To show that all the arrows of time point in the same direction we will 
consider that the master arrow of time is the cosmological one. We will 
show that the universe expansion creates a thermodynamical instability in 
the universe, in such a way that the thermodynamical arrow of time must 
necessarily point in the same direction. We will refer to the literature for the 
problem of the coincidence of the other arrows of time with the cosmological- 
master arrow. 

1.2 The problem of the unification of dynamics and 
thermodynamics. 

A particular, but very important, case of the first problem, is the prob- 
lem of the unification of the time-symmetric dynamical laws with the time- 
asymmetric thermodynamical laws. In fact, it is reasonable to think that 
thermodynamical laws could be demonstrated using the classical or quan- 
tum dynamical laws. But, it seems that this is not the case for the second 
law of thermodynamics, that says that entropy increases, in irreversible evo- 
lutions, leading the system to a state of thermodynamical equilibrium or 
maximal entropy. This problem can be state as follows: 

i.-Liouville equation is the time- symmetric evolution equation for classical 
distribution functions or quantum density matrices p. 

a. -This equation prevent the definition of any function of p : F{p) {con- 
structed only with p and mathematical elements of the Liouville-phase space) 

such that F{p)> , namely it is impossible, as a consequence of Liouville 
theorem, to define a Lyapunov variable, i.e. a growing function of p., e. g.: 
the volume or the support of a characteristic distribution function p is time 
constant, Gibbs and conditional entropies are time constant [6], etc. etc. 

Hi. -Nevertheless we actually see that the evolution leads the system to a 
thermodynamical equilibrium with a maximal entropy stationary state . 

Therefore the problem is to combine Liouville theorem with the obvious 
fact that usual physical systems have a tendency to go to a thermodynamical 
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equilibrium. The solution of the problem is based in a theorem by Mackey 
and Lasota [6] (Theorem 4.3.1 below): 

Theorem; Let S{t) be an ergodic transformation, with stationary equi- 
librium density {of the associated Frobenius-Perron operator P{t) in a 
phase space of finite p^: — measure). Then if S{t) is p^-mixing if and only if 
P{t)p is weakly convergent to p^ i.e.: 

lim(P(i)p| g)^{pAg) (1.2.1) 

t— >oo 

for all bounded measurable functions g. 

i.e. if the time evolution is phase space is S(t) and the corresponding time 
evolution of the distribution functions is p{t) = P{t)p{0), and this evolution 
is mixing, a chaotic property of evolutions that we shall define below, and if 
there is an equilibrium density such that P{t)p* — p*, then eq. (1.2.1) can 
be proved. 

But: 

lim P{t)p^p, (1.2.2) 

t— >oo 

in fact: as we shall see in many cases this limit do not even exist. Therefore 
we have a weak limit but we have not a strong limit (i.e. a limit in the norm). 

Nevertheless we never see or measure p. What we see and measure are 
mean values of physical quantities O such that: 

< O >p=< p\0 > (1.2.3) 

Thus what we actually see is that: 

^lim < O >p=< O (1.2.4) 

In fact, all the mean values of the physical quantities go to their equilibrium 
moan values if the evolution of the system is p*-mixing. So the solution of 

the problem is quite easy: 

i.-Liouville theorem is embodied in eq. (1.2.2): the system do not go- 
(strongly) toward the equilibrium states. 

a. -Tendency toward equilibrium is embodied in eq. (1.2.4): the mean 
values of all the physical quantities goes to their equilibrium values. 

Clearly these facts are not contradictory. We will call to this solution the 
non-graining solution. 
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As chaotic-mixing systems arc very frequent in the universe the problem 
is essentially solved. What it is left to be studied are the different technics 
to deal with the detail calculations. These technics try to find some logical 
modification of the theory in order to solve the missing limit (1.2.2), which, 
even if unnecessary from the mathematical point of view, it is the way physics 
used to think, (or love to think) at least up to now. In fact there are two 
technics: 

1.2.1 Cocirse-Graining. 

Let us define an arbitrary, but time independent, projector: 

P-\9){9l (5l^) = l (1-2.1.1) 

and let us define a coarse-graining density function as: 

P= Pp=\g){g\p) (1.2.1.2) 

Prom eq. (1-2.1) we have: 

hm \g){g\P{t)p) = \g){g\p,) (1.2.1.3) 

and therefore: 

limp(t)=P, (1.2.1.4) 

t— >oo 

that would be the coarsed-graining version of eq. (1.2.2) and the main equa- 
tion of the first technic (of course the same thing happens with the general 
projector 11 = XI \9i){9i\-, (fi'ilfl'j) = ^ij)- It is easy to demonstrate that (1.2.1.4 
is a limit in norm. It is also evident that eq. (1.2.1.4) can be obtained with 
a quite arbitrary state \g) and that all the philosophy, typical of the coarse- 
graining technic, namely the definition and consideration of macroscopic and 
microscopic states [8] is just an intuitive justification to give a physical mean- 
ing to the limit (1.2.1.4). But as this justification is really unnecessary, since 
the relevant and important limit is (1.2.4), the physical explanation of all 
the philosophy of the coarse-graining technic can be philosophically criti- 
cized [9]. This is the main problem with coarse-graining. It is an arbitrary 
method. It works perfectly well but it is difficult to justify based on physical- 
philosophical (metaphysical?)-arguments. 
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In fact, coarse-graining contains the miss- leading statement: we cannot 
see microscopic states (i.e. p) but we can see macroscopic states (i.e. P). This 
statement leads to the problem of finding an unique an reasonable definition 
for these macrostates. This problem is unsolved and, in our opinion, it will 
remain unsolved since \g) is essentially arbitrary. Also, if we arbitrary chose 
some definition of macrostates, we are introducing a physical element that 
really it is alien to the system itself, and therefore this definition, even if 
natural in particular examples, will be suspicious from a general point of 
view. 

The correct " no- graining" statement is: we cannot directly measure mi- 
croscopic states (i.e. p), we can only measure mean values of physical quan- 
tities or observables (among them the projector P = \g){g\ and therefore the 
arbitrarily defined macroscopic states). This statement is completely true at 
the classical and quantum levels [10] and refers to all physical observables. 
Then we can rigorously say, e.g. that the two thermodynamical variables 
< p > and < V > define the thermodynamical macrostate of a perfect gas. 
etc. 

1.2.2 Fine Graining. 

Let C is the Hilbert-Liouville space of the physical states p and = C the 
space of the hnear operator on C. We may think that not all O e C — C^ 
is a physical admissible observable. In fact, observables are measured by 
real physical devices, that very likely are free of sophisticated mathematical 
behaviors: e.g. are related with continuous and derivable functions and not 
with discontinuous non derivable functions, even if square-integrable. So it 
is reasonable that 0{x) would be ,e.g. a Schwarz function (we will precise 
this point in section 5). So let us call $ to the space of physically admissible 
observables such that: 

$ C £ = (1.2.2.1) 

If we consider the dual of $ we have a Gel'fand triplet (cf. Appendix 
4.A.): 

$ C £ = C (1.2.2.2) 

(as we shall see in section 5 if we give to the functions of $ some analyticity 
properties we can consider also the time-asymmetry problem within this 
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framework). We will work with states that belong to e.g.: normally 
belongs to this space. As it is well know we define the functional A addition 
of the functionals B and C :A — B + C as the functional A[g] — {A\g) defined 
by: 

A[g] = {A\g) = B[g] + C[g] = {B\g) + {C\g) 

for all g' G $. The same method is used to define the product of a functional 
by a number. Analogously, if we have a sequence of functionals: Ai,A2, ...the 
limit A =lim Ai is defined as the functional such that: 

i— ♦oo 

A[g] = {A\g) -lim Ai[g] =lim {Ai\g) 

I— >0O I— >oo 

for all g e 

Then, as p and can be considered as functional on $ eq. (1-2.1) reads: 
lim P(t)p^ (1.2.2.3) 

t— >oo 

and we have found a rigorous "strong" limit corresponding to eq. (1.2.1) 
[11]. Perhaps the main problem with the fine-graining technic is that it is 
usual to consider the states of \ £ as unphysical states or just effective 
states, where some characteristic of real physical states have been neglected 
(as Zeno and Khalfin effects). Nevertheless we can also say that every state 
that can be used to measure the mean values of all observables of $ is a 
physical states, and this is the case with all the states of But this point 
it is not completely clear today. 

So neither technics is completely sinless. Nevertheless as the physical real 
problem is solved by the Mackey and Lasota theorem, we can say that all this 
sins are venial sins. On the other hand both technics have some advantages: 
e. g.: 

i. -Coarse graining works just with one physical space, C. Also coarse- 
graining is unavoidable to calculate global thermodynamical variables like 
temperature or pressure, but 

ii. -The time evolution of p{t) can be computed easier using the fine- 
graining technic, since we have the vectors of space $^ that can be used 
to find new spectral expansion for the observables of the problem. Once we 
know p{t) we can compute P— Pp{t), while the direct computation of P (t) 
,using coarse-graining technics directly, can be more difficult [12]. 



8 



These lectures will be almost devoted to study the new ideas of no- 
graining and fine-graining, since the coarse graining technic is well known 



The lectures are organized as follows: 

In section two we will describe the dynamics, both classical and quantum, 
and define the notions of time-symmetry and reversibility. This section is 
based in paper [15]. 

In section three we will deal with thermodynamics and we will give dif- 
ferent definitions of entropy. This, and the two next section are based in 
paper [6] , but we have added the new mathematical and physical structures, 
recently appeared and studied. 

In section four we will introduce the classical evolution equations and we 
will study the ergodic, mixing, and, exact transformations. 

In section five will see the quantum cvohition equations. We will study 
the no-graning and fine graining ideas, both in models with discrete and 
continuous spectra, and we will consider the Friedrichs model, for pure and 
mixed states. 

In section six we will study the coarse-graining projectors and the fine- 
graining traces. We will also study the problem of time asymmetry. 

In section seven we will review the main equation of thermodynamics in 
curved space-time. 

In section eight we will consider the coordination of the arrows of time. 
This section is mostly based in references [2], [4], and, [16]. 

In section nine we will drown our conclusions. 



In this section we will review the formalism that we will use in this work 
and we will see how the notions of reversibility and time-asymmetry are 
introduced. 

2.1 Classical formalism. 

A classical system with N degrees of freedom is characterized by its Hamil- 
tonian 



(cf. [11],[12],[13]) 



2 Dynamics. 
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a function of x, the generic point in the 2N-dimensional phase space X or 
a function of the configuration variables qi and the momentum variables pi 
{i — 1,....,N). The system is solved if we compute the functions 



Qi = Qi{t), 



or 



X = x{i) 

solutions of the Hamilton equations 



(2.1.2) 




(2.1.3) 



P.(0)=P° 



(2.1.4) 



The solution of the system of differential equation s (2.1.3) is the map S{t) : 
X — > X, defined by: 



We also call St — S{t) and these St form a group. If ^4 C X is a subset of 
the phase space we can compute the image of A, namely St{A) — B. Then if 
Hl is the Lebesgue measure on X we can formulate the Liouville 

Theorem 2. 1.1. If S{t) is the map obtained solving the classical dynam- 
ical evolution and A a /Xi-measurable set of X then: 



i.e.: Classically the evolution preserve the "volumes" of phase space. 

Let us now define the notion of reversibility. Experimentally it is im- 
possible to change the direction of time. The best we can do in order to 
simulate a time inversion, is to film the motion under study and project the 
film backward. Then, if qi = qi{t) and pi = Piit) gives the real motion. 



s{t)[xm 



x{t) 



(2.1.4') 



(2.1.4") 



10 



the law of the fictitious motion obtained playing backward the film will be 
Qi = Qi{~t)j Pi = ~Pi{~i)-i where to change t by —t is simply an easy way 
to avoid to define new initial data (the final ones of the reversed motion). 
We can deduce that the time reversal operator T acts on the configuration 
variables and the momentum variables as [1],[17] 

T{qi,Pi) = {Tqi,Tpi) = {qi,-pi) (2.1.5) 

We can now consider the data (2.1.4) (that we have called "conditions 
at zero time" and not "initial conditions" in order to avoid any reference to 
time, even though we shall follow the common use in other sections) and 
compute the reversed data 

gr(0)=g,(0), 



pr(0) = -Pi{0) (2.1.6) 

With these conditions, "at zero time" we can calculate, using Eq. (2.1.3), a 
new real motion that we will call ql^'"{t), pl^'"{t). We will say that the motion 
is reversible if 

qrit)=q^i-t), 



pr{t)--Pi{-t), (2.1.7) 

namely, if the motion in the backward film agrees with a real motion with 
reversed conditions at zero time (wc see that the initial conditions of one 
motion are the final ones of the other). 

Usually H (cf. Eq. (2.1.1)) is quadratic in the pj, so that 

TH{qi,p,) = H{Tqi,Tpi) = H{qi, -pi) = H{qi,pi) (2.1.8) 

In this case we will say that the hamiltonian is time-symmetric. 
Then, if we make a T transformation (2.1.5) on Eq. (2.1.3), we find 
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_d(^^_g ^2.1.9) 

dt ' ^ ^ 

and if we now change i by — i we find again Eqs. (3.1.3) as 

d{-)Pi 



d{-t) 



d,.H (2.1.10) 



From this equation and eq. (2.1.3) a motion {q^^'" , Pi^'") with data (2.1.6) 
must satisfy (2.1.7). Therefore 

Theorem 2.1.2., A usual Hamiltonian, quadratic in the p\s, yields a 
reversible motion. 

The only condition to obtain a reversible motion is eq. (2.1.8), namely 
that the hamiltonian would be time-symmetric. 

Then reversible motion form a group. But irreversible motions do not 
form a group, since the inverse of such a motions do not even exist, because 
they are not real motions. 

We will further say that the initial conditions are time-symmetric, with 
respect to i = if Pi{0) — or: 

QiiO) = qi{0) 

Pi{0) = -Pi{0) (2.1.11) 

Then, if the motion is reversible, we will have: 

Qiit) = Qi{-t) 
Pi{t) = -Pi{-t) (2.1.12) 

We call this motion time -symmetric with respect to t = 0, since the curves 

qi{t) arc symmetric with respect to the vertical axis and the curves Pi{t) are 
symmetric with respect to the origin of the coordinate system, as in fig. 0. 
Therefore: 

Theorem 2.1.3. If the motion is reversible and the condition at i = 
is time-symmetric, the motion is time-symmetric with respect iot — Q. 

If all the motion would be time-symmetric with respect to t = it would 
be impossible to define any arrow of time at t = 0, since past and future 
would look like exactly the same from this instant of time. 
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2.2 Quantum formalism. 

The quantum wave function for the same system treated in Sect. 2.1. reads 

Hqi,t) = {qi\m) (2.2.1) 

This function belongs to a Hilbert space H — L'^. Namely if we introduce 
the inner product: 



($,^) = J <^*^^Jd^q (2.2.1') 



it is ($, $) < oo, and usually is normalized as ($, $) = 1 
and satisfies the Schrodinger equation 

i^ = Htte,0. (2.2.2) 

from which we can find the time evolution of the wave function $(gj,t) by 
imposing conditions at zero time 

$(gi,0) = $Ofe) (2.2.3) 

Then: 

^qi,t) = e-^^*$(?,,0) = um\qi) (2.3.3') 

Since we are now working in the configuration representation, in which 
the position and momentum operators are 

qi = qi, 

Pi^-idg„ (2.2.4) 

the quantum version of Eq. (2.1.5) is 

T^q,,t)^^*{qi,t) (2.2.5) 

For, if 

(Pi)^ = / $*(g)( - idq,)^q)dq, < q^ j ^*{q)qi^q)dqi 
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then 

(Pi)** = j ^{q){ - idq^)^*{q)dq = < % % >$ 

(for more details see [17]). Then, the wave function of the inverted motion 
will have as zero time data 

^reviQi: 0) = 0) = (2-2.6) 

and the motion will be reversible if 

^rev{q^,t) = ^*{q^,-t), (2.2.7) 

which is the quantum version of eqs. (2.1.7). 

If is a hamiltonian quadratic in p, it is easy to see that H is real (or 
time-symmetric) namely: 

H^H* (2.2.8) 

Then we can formulate the 

Theorem 2.2.1. If the hamiltonian is real the corresponding evolution 
is reversible. 

Proof. : 

From eqs. (2.2.6) and (2.2.8) we can obtain (2.2.7), since 

^rev{t) = e-^^*$,e.(0) = e-^^*$*(0) = (e'«*$(0))* = $*(-t), (2.2.9) 

where we have omitted the variables qi. Then, as in the classical usual 
Hamiltonian yields a reversible motion. □ 

We can also show directly that Eq. (2.2.2) is t-invariant, but we prefer 
the proof above because the role played by the condition at zero time can be 
seen explicitly. 

As in the classical, case reversible motion form a group, since u~^{t) = 
u{—t) is a real motion, which is not the case for irreversible motion where 
this motion is not a real one. If u{ti)u{t2) = u{ti + ^2) for ti, t2> only, we 
will say that these motions form a semigroup. This is the case for irreversible 
motions. 

Let us now repeat all this formalism, that we have so far introduced in 
the configuration representation, in an abstract way. The state of the system 
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is defined by the ket that belongs to the Hilbert space of states Ti, 

and satisfies the Schrodinger equation 

ij^\m)-H\m). (2.2.10) 

The inner product is symboUzed as < >= and the normahzation 

is < $1$ >= 1. 

> can be found solving eq. (2.2.10) with the condition at zero time 
|$(0)) = (2.2.11) 

Namely: 

|$(t) >= e-^-^*|$(0) >= M(t)|$(0) > (2.2.11') 

Then, the T transformation can be defined as [17], [18] 

T\^{t)) = K\^{t)) = \^*{t)), (2.2.12) 

which means that we must conjugate the wave function in the configuration 
representation and then go to the generic representation. K is known as the 
Wigner operator. More precisely, let {^{qi, t)) — \^{qi, t)) be the coordinates 
of the state vector in the configuration representation (wave function) and 
|$(, i) > the coordinates of the same vector in a generic representation; then 

|$(t) >=t/|$(g„t)), UU^ = 1 (2.2.13) 

Let Kq be the conjugation operator in the configuration representation 

i^0|$(g^,t)) = |$*(g^,t)); (2.2.14) 

then 

K\^{t)) = KU\^{q,,t)) = UKo\^{qi,t)) = UKoU^\^{t)) (2.2.15) 

Namely, if Kq is the conjugation in the configuration representation, the 
Wigner operator X in a generic representation reads 

K^UKoU^ (2.2.16) 
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It is easy to show that in the configuration representation we have Kq 
has the following properties: 

a) Kq is an antilinear, antiunitary operator, namely [17]: 

al) Koia\l) + P\2)) = a*Ko\l) + P*Ko\2); 

a2) if |2) = Ko\2), (1| = {1\kI and A = KqAKI then (i|i|2) = 
(11^12)*; 

a3) ((l|A'o)|2) = (l|(i^o|2)r, 
i.e., parentheses cannot be omitted. 

b) Kq — 1 (at least for spin zero fields [17]); 
c) 

KopiKl = -pi, 

d) 

KqcKI^c* if ceC. 

Therefore, KqKq = 1 and KqiKq = — i. 

From (2.2.16) it is also easy to show that K has the same properties. 

As an exercise we can repeat formulae (2.2.6) to (2.2.9) in a generic rep- 
resentation. The time reversal is given by Eq. (2.2.12). The reversed initial 
condition is 

mO)rev) =KmO)), (2.2.17) 

and the condition of reversible motion reads 

mt)rev) ^ Km-t)) (2.2.18) 

We will say that H is real if 

H^KHK\ (2.2.19) 

and usually H is endowed with this property, because Eq. (2.2.8) is satisfied 
in the configuration basis. Then, from Eqs. (2.2.17) and (2.2.19) we can 
deduce again theorem 2.2.1, now in a generic coordinate system: 

mUv) = e-'^*|$(0),e.) = e-^*i^|cl>(0)) = 
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= K{K^e-'"^K\^{0))) = Ke'"^\^{0)) = K\^{-t)) (2.2.20) 

Then, as in the classical usual real Hamiltonian yields a reversible 

motion. 

In general, we will call a kct |1) (bra (1|) real if 

m = \^), or 

{1\K^ = (1|, (2.2.21) 

and an operator A real if 

KAK^ = A (2.2.22) 

From (2.2.19) we see that a usual Hamiltonian is a real operator. 
A basis {|i)} will be a real basis if all its kets are real: 

K\i) = \i) (2.2.23) 

In a real basis, K is just the conjugation of the coordinates of the vectors, 
or of the coordinates of the operators: 



K\c|>) = KJ2c^\^) = Y: 



C, 1} 



KAK^ = K{Y.c,^\i){j\)K^ = E4N)01 (2-2.24) 

Therefore, the configuration basis is {\x >} real. 

We will say that the conditions at i = 0, are time-symmetric if: 

|$(0) >= K|$(0) > (2.2.25) 

namely |$(0) > is real. Then, if the evolution is reversible we have: 

\^{t) >= K\^{-t) > (2.2.26) 

and we will say that the evolution is time-symmetric with respect to t = 0. 
So we have the 

Theorem 2.2.2. If the evolution is reversible and the initial condition 
is time-symmetric, the evolution is time-symmetric. 

Then we can repeat what we have said in the classical case. If all the 
quantum evolutions would be time-symmetric with respect to i = it would 
be impossible to define a quantum arrow of time at t = 0. 
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2.3 Statistical formalism. 



We shall treat simultaneously the classical and quantum cases in order to 
establish an analogy or unified formalism that we shall use below. Never- 
theless, it must be clear that there is a great difference between the classical 
and quantum cases. 

We will call the classical distribution function or density (resp., the quan- 
tum density matrix) a function (resp., matrix) endowed with the following 
properties: 

p{Qi:Pi) > 0; ..or...p{x) > 

II P 11= / p{Qi:Pi)dqidpi = / p{x)dx = 1, (2.3.1) 

where X is the phase space. Distribution functions p belong to a Hilbert 

space called the classical Liouvillc space. 

(resp., in the quantum mechanical formalism: 

P = p^; 



tr{p) = 1; 

Paa > 0) (2.3.2) 

Density matrices p belong to a space JC — H x H called the quantum Liouville 

space.) 

p satisfies the Liouville equation 

idtp = Lp, (2.3.3) 

where 

L = i{H,..}pB (2.3.5) 

(resp., 

L=[H,..] = H xl-lx H (2.3.5) 

(cf. eq. (2. A. 24) for the definition of x)). 
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Therefore the time evolution, in both classical and quantum cases, is: 

= e-^V(O) = C/Wp(0) (2.3.5') 

The T transformation of a density function is 

Tp(qi,Pi) = p'iQhPi) = p(qi, -Pi) (2.3.6) 

(resp., the T transformation of a density matrix is 

Tp^p' ^ KpK^ = /Cp (2.3.6'). 

where IC = K x K^). 

Prom Eq. (2.3.3), if the Hamiltonian is a usual time-symmetric one, we 
have classically 

TLp{qi,Pi) = Ti{H,p}pB = 
Ti dq^Hdp.p - dp.Hdq.p = 

i 

iY.d,m_pTp-d_p^Hd,Tp^ 

i 

-i{H,Tp}pB = -LTp; (2.3.7) 

therefore, if we T-transform classically Eq. (2.3.2), we obtain 

idtTp= -LTp (2.3.8) 

(resp., if we T-transform the quantum Liouville equation (2.3.2), we obtain, 
if the Hamiltonian is a real usual one, 

KiK'^dtKpK'^ ^{K X K'^)L{K'^ x K)KpK^; (2.3.8') 

but KiK^ = -i, and {K x K^)L{K^ x K) = KHK^ x 1 - 1 x KHK^ = L, 
so 
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(2.3.8") 



i.e., the same equation as the classical one (2.3.8)). 

In both cases a minus sign appears. In the reverted solution we must 
change t by —t, namely: 

T{t) =t' ^ -t (2.3.9) 

So we have proved the 

Theorem 2. 3.1. The Liouville equation remains invariant under T trans- 
formations for a usual time-symmetric Hamiltonian. 

Thus we have shown the complete isomorphism of the classical and quan- 
tum formalisms. From now on we will mainly use the quantum formalism, 
since it is the one that is better known by physicists. Let us therefore review 
the main properties of the usual Hamiltonian, in a real basis to simplify the 
treatment. Prom the equations of Sect. 2.2. we have 

H^H^, (2.3.10) 

namely, the Hamiltonian is 

I- self- adjoint, because it is an observable; 

II- real, because for the usual Hamiltonian the motion is reversible; 

III- as a consequence, it is also symmetric. 

p belongs to a set that, endowed with the inner product (2.A.1) of the ap- 
pendix, becomes the Liouville-Hilbert space C From Eqs. (2. A. 29) and 
(2.3.10)) we can prove that the Liouvillian has the following properties in 
real basis: 

L = L\ 

L = L\ 
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L^-L^ (2.3.11) 



Then, from (2.3.11) and (2. A. 18) we have 

(iL) = {iiy (2.3.10) 

This property is important since, from it, wc can deduce that the matrix 
p remains Hermitian under the evolution satisfying to the LiouviUe equation 
(2.3.2). In fact, it follows from Eq. (2. A. 18) that a product of self-associated 
commuting operators is also self-associated. Then, {iL) = {iLY implies 
^^-iLty = (,-iLt ^^Q-) ^ p(^Qy implies e-^^V(O) = [e-'^^p{0)]\ namely. 

Pit) = Pity. 

Finally, let we prove once more that if the Liouvillian is real the evolution 
is reversible. Based on Eqs. (2.1.7) and (2.2.7) we define a reversible motion, 
in a real basis, as 

Prev{t)=P*i-t), (2.3.13) 

where Prev{t) is the motion with reversed condition at zero time: 

Pre.(0)=p*(0) (2.3.14) 

Now we can prove the 

Theorem2.3.2. If the liouvillian is real the evolution is reversible. 
Proof. : 

Then, with the same reasoning as for eqs. (2.2.9) we have: 

Previt) = e-^Ve.(0) = e-^V(O) = [e^^V(O)]* = P*{-t), (2.3.15) 

which shows that a motion with a real Liouvillian is reversible □ 
.In a generic basis eqs. (2,3,11) tO (2.3.15) read as follows: 
The liouvillian is real or time symmetric if:. 

KLlO = L, (2.3.16) 
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The evolution is time-symmetric if: 

Previt) = Kp{-t)K'^ = ]Cp{-t) (2.3.17) 

The conditions at time t = are time-symmetric if: 

PreM = Kp(0)Xt = /Cp(0) (2.3.18) 

A real liouvillian and time-symmetric conditions at t = yields a time- 
symmetric evolution since: 

Prev{t) = e-*^Vre.(0) = e-''^' K p{Q)K^ = [e*^V(0)]i^^ = Kp{-t)K^ (2.3.19) 

The condition at i = will be called time-symmetric if: 

p(0) = Kp{Q)K^ = ICp{0) (2.3.20) 

then if the evolution is irreversible we have: 

p{-t) = ICp{t) (2.3.21) 

and we will say that the whole evolution is time symmetric and we can repeat 
what we have said in the previous cases. So we have the 

Theorem 2.3.3. If the evolution is reversible and the condition at t = 
is time-symmetric the evolution is time-symmetric with respect to t — 0. 

Proof.: 

If the liouvillian satisfies eq. (2.3.16) and condition at t — satisfies eq. 
(2.3.20) all the evolution is time-symmetric, since: 

}Cp(t) = ;C(e-*^V(0)) = e^^'-^^'JCpiO) = e^^V(O) = p{-t) (2.3.22) 

Therefore the motion is time-symmetric if L is real and the condition at time 
t = is time symmetric. □ 

2.4 Appendix 2 A. Mathematical theory of superspace 
and superoperators [19]. 

Let us make a small mathematical interlude, to define the notions of super- 
space and superoperators. 
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2.4.1 The quantum case 

Let us consider a Hilbert space Ti and the space £ = x 7i of matrices on 
7i, i.e., the Liouville-Hilbert space. Matrices will be symbolized by greek 
lower case letters a, P,....,p, with coordinate a^j, Pij,....,pij. We will call the 
linear space of matrices the superspace L and the matrices supervectors. Let 
us define an inner product in the superspace JC: 

a-P= {a\P) = tr{a^P) = J^a*jPij (2.A1) 

ij 

Using this inner product L becomes a Hilbert space. 
The norm of a supervector is thus 

||q;|| ^a-a^J2M'^ (2.A2) 

ij 

Let us consider the hnear operators in superspace, that we shall call 

superoperators, and that we shall represent by capital Latin letters A, B, 
L, with coordinates Aij^^h Bij^ki- Superoperators act on matrices as 

Aa = P, 

aA = (3 (2.^.3) 

We will use, for these two equations, the following rule for indices 

^^ij,kl<^kl = l3ij, 
kl 

^alA,,, = pJ, (2.A4) 

kl 

In the first equation we have used the usual multiplication "row by col- 
umn" and a and /3 are considered as column vectors. In the second one we 
have transposed a and /3 since, in these case, they are considered as row 
vectors. 
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Since the supcropcrators have four indices we can define more operations 
defining transposed and adjoints than for ordinary two-indices matrices. So, 
we define, for a superoperator A, its 

a- Transposed as the superoperator such that 

Aa = aA^ {2.A.5) 

for aU a & C. 
Then, 

Hk,ji 

SO 

'kj,kl = ^lk,ji 

Of course. 



Aj,ki(^ki = oikiAjuji, (2.^1.6) 



Aij,ki = ALi (2-A7) 



{A^f = A, 

{A,A2f = AlA\ (2.A8) 

and A is symmetric (antisymmetric) if 

A = A^, {A = -A^) {2.A.9) 

b- Adjoint A^ as the superoperator such that 

Aa = {a^A^y (2.A10) 

for all q; e £. Then, 

A^JMakl = = MA,^jT, (2.^.11) 

so 

^ij,kl = ^kl,ij (2.74.12) 

Of course, 

(A)t = A, 
(A.A^)^ = {AIa\), (2.A13) 
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and A is Hermitian (anti-Hermitian) if 

A^A\ (A^-A^) (2.^.14) 

c- Associated A'' as the superoperator such that 
Aa = (AV)^ (2.A15) 

for all q; e Then, 

Aij,kiaki = iA^iaiy = {A^M^^-^.y = {A^ ai^), (2.^.16) 

so 

Aij,kl = '^'ji,lk (2.74.17) 

Of course, 

= A, 

{A1A2)'' = A^A'^, (2.A18) 

and an operator is adjoint- symmetric (or self-associated) if 

A = A (2.^.19) 

An adjoint-symmetric operator acting on a Hermitian matrix gives an- 
other Hermitian matrix. For, if 

a^a\ A^A\ (2.A20) 

then from Eq. (2. A. 15) we have 

Aa = (Aa)^ (2.A21) 

Putting all together, we have 

AjM = = (Al-.r = (A« (2.^.22) 

and therefore 
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= (2.A23) 

Let us now define a superoperator as a product of two operators, A — 
Q! X in the following way: 

A-f^a-f(3, V7, (2.A24) 

or, equivalently, 

kl kl 

that is, 

Aij,ki = aikPij (2.yl.25) 

Then, 

^IkiAik.ji = ^1kiOiijl3ik = ^PiklkiOiij (2.^.26) 

kl kl kX 

and, from (2. A. 4), 

-fA = (2.A27) 

Therefore, we have from (2. A. 24) and (2. A. 27), 

{a X /3)7 = q;7/3, 

7(a X /?) = (2.A.28) 

The choice of the index position, in eq. (2. A. 4), was made in order to 
obtain these simple multiplication rules. It is easy to prove that 

{a X (3Y = /3 X ct; 

(a X /3)f = at X /?t; 
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{a X = /3t X (2.A29) 

The product x can be used to define the time inversion of matrices, since 
a time inverted matrix is (Eq. (2.2.22))): 

Tp = KpK'' = (X X K'')p = JCp (2.A.30) 

From this equation we can deduce the time inversion rule of superoperators, 
namely, 

TA^{Kx K^)A{K X K^)^ = /CA/C^; (2.A31) 

since (cf. Eq. (2.A.29)) 

{K X ^K^ xK, 
we have the alternative expression 

TA^{K X K^)A{K^ X K) (2.A32) 

We can also compute {a x P)^^ x 6): 

J2{01 X P)ij,kl{l X S)kl,nm = 
kl 



(^ikPljlkn^ml — X] ^iklkn^mlPlj — 
kl kl 



(q;7 X 5(3)ij,nm, (2.^.33) 

namely, 

{a X /3)(7 X 5) = {a-f x 5(3) (2.A34) 
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2.4.2 Classical case. 



As wc have seen the quantum Liouville space is transformed in a Hilbert 
space by the inner product (2.A.1). In the same way it is convenient to define 
an inner product in the classical Liouville space L ,namely: 

(p|(7) = / p'{x)a{x)dx (2.A35) 

Using this inner product and Wigncr functions (App. 6. A) the classical C 
becomes also a Hilbert space, and the classical analogs of the quantum 
equation of the previous subsection can be found. Also we can use the Wigner 
function integral of appendix 6. A to make the analogy explicit. 



3 Thermodynamics. 

3.1 Classification of the different types of second laws. 

The first law of thermodynamics is just the conservation of the energy. There 
is no confiict between dynamics and thermodynamics about this law. The 
problem is to derive the second law of thermodynamics based in dynamical 
considerations. The second law is expressed in many forms by the different 
authors, therefore we will begin our research by a classification of these forms. 
Let S{t) denote the thermodynamical entropy of a closed system: 
i.-We will call a first-order second law to the statement: 



S{t) > S{t') (3.L1) 

if t > t', thus according to this form the entropy cannot decrease. 

ii.-An stronger assertion would be a second-order second law: Eq. (3.1.1) 
is satisfied and also: 



lim S(t) = (3.1.2) 

in this case we assert that the system entropy converges to a steady-state 
value S^, which may not be unique, e.g. it can be the entropy of a metastable 
state. Different preparations of the system could yield different final metastable 
states. 
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iii.-The final and stronger form, or third-order form, of the second law is: 
Eqs. (3.1.1) and (3.1.2) are satisfy but also the limit (3.1.2) is unique. In this 
case the entropy of the systems evolve to a unique maximum, irrespective of 
the way the system was prepared. 

We will find these different forms of the second law below. 



3.2 Dynamics and densities. 

We will consider more generic systems than the ones of section 2, in order to 
be as general as possible. This is not done just by the sake of mathematical 
generality but because we will be forced to consider such a systems below, 
if we want to solve our problems. So let us consider a system operating in 
a phase space X, with a evolution law St more general than (2.1.3), i. e. 
a mapping St : X ^ X, that change the point a; of X as t changes. X- 
may have finite dimension d or infinite dimension, and t can be discrete or 
continuous. We will consider only "autonomous" processes: i.e. such that 
St{St'{x)) = St+t'{^), So{x) = X. Thus the mapping S can form either 
a group of transformation when t' e i?, (or Z) ( e. g. the evolutions 
with time. symmetric hamiltonian or liouvillian of sec. 2) or a semigroup if 
t,t' G -R"''(or A^). In the two last cases {R'^,N) an equation like (2.1.7) do 
not exists and the evolution is necessarily irreversible. 

For every point xq the successive point St{xo) are a system trajectory. To 
study an infinite number of initial point or an infinite number of trajectories 
we introduce the density functions p{x) e L^{X) , namely: 



and such that: 
where: 



/ \p{x)\dx <oo, (3.2.1) 

pix)>0,\\ p{x) Ihl, (3.2.2) 

II p[x) 11= / \p{x)\dx (3.2.1) 



is the L^-norm of p. We use to postulate that a thermodynamics system is 
a system that has , at a given time, states distributed throughout the phase 
space X and the distribution of these states is characterized by the density 
function p{x). 
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We will call the p-measured lJ>p{A) of the set ^4 C X to: 
fXp{A) = / p{x)dx (3.2.4) 

J A 

The Lebesgue (non-normalized usual) measure of a set 74-will be denoted 
Hl (A). The uniform density will be 

Pl{x) = (3.2.5) 

and therefore Lebesgue normahzed measure is iJip^{X) = ij,l{X) — l.We 

always write /il (dx) — dx. 

Finally X can be either Gibbs phase space F or Boltzmann phase space 

3.3 Gibbs entropy. 

This entropy is defined as: 

H(p)^^[ p{x)logp(x)dx (3.3.1) 

Jx 

It is and additive quantity, namely the Gibbs-entropy of a system form by 
two subsystem is the sum of the two corresponding entropies. Then it is 
called an extensive quantity. Gibbs-entropy can be written as: 

H{f) = I v{p{x))dx (3.3.2) 

where the r/(p) function is defined as: 

r]{p) — — p log p.. /or.. p > Q..and..r]{Q) — 0.... (3.3.3) 
and it is endowed with the property: 

r}{p)<{p-a)r^'{a)+rj{a) (3.3.4) 

combining these last two formulae we can prove the Gibbs inequality: 

p — plogp < o" — plogo", ..for..p, a > (3.3.5) 

If p and (7 are two normalized density function, integrating the last equation 
we have: 

— / p{x) log p{x)dx < — p{x)\oga{x)dx..{3.3.6) 
Only when p = a does the equality hold in eqs. (3.3.4), (3.3.5), and (3.3.6) 
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3.4 Microcanonical and canonical ensembles. 



Let us consider an space X with a finite Lebesgue measure: (X) < 
oo.Then the only density that will make Gibbs-entropy maximal is the uni- 
form density of eq. (3.2.5). Precisely: 

Theorem 3.4.1.When /^l (X) < oo the density that maximized the 
Gibbs-entropy is the uniform density, Pl(x) (cf. eq. (3.2.5). For any other 
density p 7^ PlH{p) < H{pl). 

Proof.: Chose an arbitrary density p ,thus from eq. (3.3.6) we have: 

H{p)<-f p{x)\oga{x)dx (3.4.1) 

However, if a{x) = 1/ A*l {X) the integrated Gibbs inequality (3.3.6) gives: 
H{p)< -log (3.4.2) 

./^L (X). 

since p is normalized to one. The equality holds if p = , but the entropy 
corresponding to is: 

//(pi) = -log (3.4.3) 

therefore i?(p) < H^Pl) for any density p and i?(p) < H{pl) for p ^ p^. 

Clearly if X is normalized so /^l (X) = 1, then H[p) < 0. □ 

The uniform density it is also called the density of a microcanonical en- 
semble, and, as we can see, to define it we do not need to use any particular 
property of the thermodynamical system under consideration. 
Another, even more interesting theorem is the following: 
Theorem 3.4.2.Assume that a non-negative measurable function a{x) 
is given as well as an average or expectation mean value < a >p of that 
function over the entire X, weighted by the density p: 



< 



a >p= / a{x)p{x)dx (3.4.4) 



Then the maximum of the Gibbs-entropy H[p) subject to the constraint 
< a >p— const, occurs for the density: 



p*{x) 



.(3.4.5) 
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where: 

Z = I e-^"(^)cia; (3.4.6) 

Jx 

and i^-is implicitly defined by the normalization condition: 

< a >p= Z-^ f a(x)e-'^°^(^)dx (3.4.7) 

Proof: The proof again uses the integrated Gibbs equality (3.3.6) so: 

H{p) < - Jx p{x) log p^{x)dx = - Jx p{x)[- log Z - iya{x)]dx ^ (3 4 8) 
= log Z + u Jx p{x)a{x)dx = log Z + u < a >p 

However, it is equally easy to show that: 

H{p,) = logZ + u <a>p (3.4.9) 

and therefore H{p) < H{p^), with the equality holding if and only if p = p*, 

□ 

If a{x) is the energy of the system is density of the Gibbs canonical 
ensemble at temperature T — . (With many constraints < ctj >p we 
would define the density of a grand canonical ensemble). 

We use to postulate also that there is a one-to-one correspondence be- 
tween thermodynamical equilibrium states and the states of maximum en- 
tropy. Then, from the last theorems it would be natural to postulate also 
that the thermodynamical entropy S coincide with Gibbs entropy H{p). In 
fact with this postulate we can obtain usual equilibrium thermodynamical. 
But, as we shall see below, this identification is not what we need to base 
a non-equilibrium thermodynamics, since Gibbs-entropy has not the wright 
properties, in this case. 



3.5 Reversible and irreversible systems. 

In section 2 the properties of the hamiltonian force the motion to be either 
reversible or irreversible.. But in this section we are studying more general 
motions so we are force to repeat these definition for this more general cases. 
Nevertheless, in order to prove some theorems, the motions cannot be com- 
pletely general so we will restrict ourselves to motion produced by Markov 
operators. 
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Any linear operator: Pt : ^ such that: 



{a)..Ptp > 0, 



(b): II 



p 



(3.5.1) 



for all t e R, p > 0, p e is a Markov operator, i.e. an operator that 
acting on a density gives a density. Markov operator have a number of useful 
properties. The most important is that if p e and it is not restricted to 
p > 0, then: 



for all t,t' E R (or Z ), namely reversible Markov operators form a group. 
Evolution operator U{t) = e*'^* of eq. (2.3.5') is an example of reversible 
Markov operator. It is so because it is generated by a time-symmetric or real 
liouvillian L. 

However, if in the last definition we substitute R and Z by R^ and we 
have the definition of an irreversible Markov operator. Irreversible Markov 
operators form a semigroup. 

Gibbs entropy cannot be use in non-equilibrium theory since it may de- 
crease under the action of some Markov operators (cf. [6]), therefore we 
cannot use this entropy to formulate a second law of thermodynamic, even 
in the first-order form. Nevertheless Gibbs-entropy is completely successful 
in equilibrium situations, so the entropy we will chose, for non-equilibrium 
situations, must coincide with Gibbs entropy at equilibrium. 

3.6 Conditional entropy. 

If p and a are two densities such that suppp Csuppa, then the conditional 
entropy of density p, with respect to density a, is:(3.6.2) 



PtP \\<\\P 



(3.5.2) 



which is known as contractive property. 

A Markov operator is reversible (or time-symmetric) if: 



{a)..Pop = p, ..{b)..Pt{Pt'p) = Pt+t'P (3.5.3) 




(3.6.1) 
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The conditional entropy is always definite, i. e.: it is finite or equal to -oo. 
As it is evident the conditional entropy measured the deviation of density p 
from density a. 

Conditional entropy has two very important properties: 

i. -Since p and a are both densities the integrated Gibbs inequality (3.3.6) 
implies that Hc{p\o') < . It is only when p = a that the cqiiality hold. 

ii. -If Pl is the constant density of the microcanonical ensemble throughout 

the phase spa,ce X then Hc{p\pl) = H{p)— log pi {X). Therefore conditional 
entropy is a generalization of Gibbs entropy in this case. 

As Hc{p\p*) — when p — p^ it is reasonable to postulate that: 

S - = Hc{p\p.) (3.6.2) 

e.g. when p^ is the density of the canonical ensemble. We will see that this 
definition is completely satisfactory and that using equation (3.6.2) we can 
formulate the second law of thermodynamics in its second and third-order 
forms. The first result, along these hnes, is a weak, first-order form, of the law 
of thermodynamics, namely that the conditional entropy is never decreasing, 
as it is proved by the 

Theorem 3.6.1. [20] Let Pt be a Markov operator. Then: 

HciPtp\Pt<j) > Hc{p\a) (3.6.3) 

for all densities p and a. 

A second result is the following: if a = p^, is stationary, namely Ptp^ = p^ 
, then: 

Hc{Ptp\p*) > Hc{p\p.) (3.6.4) 

Thus this conditional entropy is always a non decreasing function bounded 
above and i?max = Hc{pJ\p*) — 0. Therefore this conditional entropy con- 
verge as t ^ oo , though more information about the evolution is required 
to find the limiting value. Furthermore if the stationary density is uniform, 
namely the one of the microcanonical ensemble we have: 

H{Ptp) > H{p) (3.6.5) 

for all non- negative p. Now i/max = ~log[l/ I^l (-^)] and as, in the general 
case we have convergency when t ^ oo. 
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Therefore eq. (3.6.2) seems a reasonable assumption. But when the 
Markov operator is reversible all these nice inequalities become equalities 
and the problem of the thermodynamical entropy reappears. In fact: 

Theorem 3.6.2. If Pt is a reversible Markov operator, then the condi- 
tional entropy is absolutely constant for all times t and equal to the value 
determinated by the choice of the initial densities p and a. That is: 

Hc{Ptp\Pt<j) ^ Hcip\a) (3.6.6) 

for all t. 

Proof.: Since Ft is reversible, by the previous theorem it follows that: 

Hc{Pt+t'P\Pt+t'(r) = Hc{PtPt'P\PtPt'o) 
> Hc{Ptp\Pt(j) > Hc(p\a)...{3.6.7) 

for all t, t' since Pt is reversible. So let us chose t' = —t, then for all times we 
have: 

Hc{p\a) > Hc{Ptp\Pt<j) > Hc{p\a) (3.6.8) 

and therefore: 

Hc{Ptp\Pta) = Hc{p\a) 

for all t. □ 

So in this case the conditional entropy is, for ever, fixed and determined 
by the method of preparation of the system. So we have gain nothing if the 
Markov operator is reversible. 

3.7 Appendix 3. A : The physical interpretation of non- 
equilibrium Gibbs-entropy. 

Gibbs-entropy is completely successful in equilibrium cases. Therefore, even 
if not the correct definition of entropy in non-equilibrium cases, it must have 
some physical meaning in these last cases. In fact, to create an unstable 
state, with a decreasing of entropy, it is necessary to provide some energy to 
the system, on the other hand, if a systems evolves from an unstable state 
toward equilibrium it release energy and the entropy grows. Therefore there 
we must find some relation like 

AS ^ -CAE + C , C > 0, (3.A1) 



3.6.9 
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This is the case since from eq. (3.6.2) we can obtain: 

AS = AHc = - / p2log— + / pi^dx (3.A2) 

J p* J p* 

where pi and p2 are the initial and final distribution functions. Now let 
us suppose that the final state is a canonical ensemble equilibrium state at 
temperature T (cf.eq.(3.4.7). Then we obtain: 

AS^AH-^J{P2- pi)ujdx ^AH- ^ (3.A3) 

so we have obtained the equation we were looking for: in fact, there is a 
relation like (3.A.1) where the coefficient is C = and C = AH.. We 
have also obtained the interpretation of the variation of the Gibbs entropy in 
irreversible evolutions: it is the difference between A^* and the linear term 
T~^AE, in such a way that if AH — 0, all the energy is used to produce 
entropy or all the entropy is used to produce energy. Therefore the difference 
of Gibbs entropy measure, some how, the efficient of the system to produce 
entropy or energy. On the other hand, if the system is isolated, so AE = 0, 
and if the equilibrium state corresponds to the one of a canonical ensemble, 
Gibbs entropy coincide with thermodynamic and conditional entropy. 

4 The Classical Evolution. 

In this section we will study "classical evolutions", in the sense that these 
evolution are not quantum ones. Nevertheless the evolutions will be as gen- 
eral as the one of the previous section: i. e., not necessarily those of section 
2. 

4.1 The Frobenius- Perron operator. 

A transformation St is called a measurable transformation if p^^^S^^ (A)) is 
well defined for all subsets A C X, where S^^{A) ~ B is the counterimage of 
Anamely: St{B) — A. Let us remark that even if a unique S^^{x) may not 
exist (as ion the case of irreversible evolutions) the counterimage do exists 
since it is the set of all the points x & B that will go to A under the action 
oiSt. 
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The transformation is non-singular if yU*(S'f ^{A)) = 



fj,{A) = 0. 



If St is a non-singular transformation, then the unique operator: Pt : 
— > L} defined by: 



is called the Probenius- Perron operator corresponding to St- For each St the 
Probenius- Perron operator is unique. If p > then Ptp > 0. As S^^{X) — X 
then II Ptp 11 = 11 p II and these operators are Markov operators. Operator 
U{t) of eq. (2.3.5') is a Frobenius-Perron operator. 
St is p-measure preserving if: 



for all sets A. Measure- preserving transformation are necessarily non 
singular, wc will also say that the measure pp is invariant under the trans- 
formation. Liouville theorem shows that transformation U {t) of eq. (2.3.5') 
is Lebesguc-mcasurc preserving. 

We will call a state p steady if PfP = p, for all t. We will call it also a 
state of "thermodynamical equilibrium", and we will symbolize it by p*. 

The relation between invariant measures and Probenius-Perron operator 
is stated by the 

Theorem 4. 1.1.- [6]. Let St be a non singular transformation and Pt 
its Frobenius-Perron operator. Then there exist an state of thermodynamic 
equilibrium whose density is an stationary state of Pt if and only if the 
measure p*: 



is invariant. 

Therefore transformation U (t) that preserve Lebesgue measure, has neces- 
sarily a equilibrium steady state, e. g. the uniform state of the microcanonical- 
ensemble. But theorem 4. 1.1. says nothing about the uniqueness of the equi- 
librium state. We shall discuss this problem in the next section. 

A point X & A a X it is called a recurrent point if there is some time 
t > such that St{x) e A. An important result is Poincare recurrence 

Theorem 4. 1.2. -Let St be a transformation with an invariant measure 
p* operating in a finite space X, p*(X) < oo, and let A be a subset of X 
with positive p*-measure. Then there exists a point x in A that is recurrent. 




(4.1.1) 



f^,iSi\A)) = ppiA) 




(4.1.2) 
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Proof.: Assume the contrary, i. e. that there are no recurrent point in 
A. This then imphes that Sj^^{A) fl A = for all times t > 0, and thus 
that S^^{A) n S^^{A) — for all positive times t ^ t' . However, since St is 
measuring preserving, this implies that iJi^{Sl~^{A)) = fi^{S^^{A)) and this, 
coupled with the pairwise disjoint nature of the sets S^^{A) and S^^{A) leads 
to: 



j:^^M)=j:^^*is;\A)) 

t=0 t=0 



[jSi\A) 
.t=o 



< /i,(X) < oo (4.1.3) 



The only way in which this inequality can be satisfied is for ii^{A) to be zero, 
which is a contradiction. Thus we conclude that A contains recurrent points. 

□ 

Therefore, in ordinary mechanical motion of finite systems almost any 
point is recurrent, since the sets of non recurrent point have measure zero. 
This fact seems to prevent the existence of irreversible evolutions, namely it 
is impossible to reach a final equilibrium state, since the system will came 
back as closed to its initial condition as wc wish, if we wait enough. The 
time wc must wait it is called the Poincare recurrence time. There are two 
ways to avoid this problem: 

i. -The practical way is to compute the recurrence time. It turns out that 
in usual system (say with a number of molecules of the order of Avogadro 
number) the time if much bigger than the age of the universe, so the returning 
to the initial conditions is practically unobservable. 

ii. -The theoretical way is to consider that irreversibility is not a notion 
of classical mechanics, but a notion that can only be defined in statistical 
mechanics, where we deals with statistical ensembles of identical system. 
Then the recurrent time of the ensemble, namely the time such that we 
reobtain the initial condition in each one of the infinite identical systems is, 
of course, infinite and the problem is theoretically solved. 

In the following subsections we shall study some properties of dynamical 
system ordered by their increasing chaotic behavior. 



4.2 Ergodicity. 

It would be interesting to know if the equilibrium state of theorem 4.1.1 is 
unique or not. To answer this question we must introduce some new concepts. 
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i. -A set A such that Sf^{A) = A is called an invariant set. 

ii. -Any invariant set A such that = or \ A) = is called 
trivial. 

iii. -A non-singular transformation St is called p*-ergodic if every invariant 
set A is a trivial subset of the phase space X. 1. e., either = or 
/i*(X\ A) = 0. This means that, if we consider a generic non-singular subset 
A, the time evolved counter-image of this subset, Sf^^{A), will wonder around 
all X since A cannot be invariant. 

iv. -If is the uniform density of the microcanonical ensemble we will say 
that St is uniformly ergodic. 

The motion within almost all tori of integrable classical mechanical sys- 
tems is ergodic [21],[22].Ergodicity is therefore a very usual property of the 
mechanical systems of section 2. 

The connection on the uniqueness of the equilibrium state and the prop- 
erties of the operators is stated in the following 

Theorem 4.2.1.[6].-Let St be a non-singular transformation and Pt the 
corresponding Frobenius-Perron operator. St is p*-ergodic if and only if Pt 
has a unique state of thermodynamic equilibrium with associated stationary 
density p*, namely a density such that Ptp* = p*. 

Hence ergodicity is the necessary and sufficient condition for the unique- 
ness of thermodynamical equilibrium, and allows us to formulate a third- 
order form second law. But this is of course half the picture, because we 
must also understands why the system evolves to this equilibrium state. 

Let us sate an important 

Theorem 4.2.2. [6]. Let St be a non singular transformation and Pt the 
corresponding Frobcnius-Perron operator with stationary density p* > for 
all points in phase space X. Then St is p*-ergodic if and only if {Ptp} is 
Cesaro convergent to p* for all densities p. i.e., if and only if 



1 



t-i 




k=0 



(4.2.1) 



in the discrete time case, or if and only if: 
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in the continuous time case, for all bounded measurable functions a and 
where (p, a) = Jx p{x)a{x)ii{x)dx (in this case n{x) is an arbitrary measure) 
is a generalization of inner product (2.A.2.1) 

4.3 Mixing. 

This will be the main property of dynamical system that we will study, it 
serve to guarantee the approach of the system to an equilibrium state. 

Let St be a p^-mcasurc preserving transformation operating in a normal- 
ized space X {ii^{X) — 1). Then St is called p*-mixing if: 

lim ii,{Ar]S-\B)) = iJ,{A)^i,(B) (4.3.1) 

t — >oo 

for all sets A and B. If is the uniform density of the microcanonical 
ensemble then we will say that St is uniformly mixing. 

Some tori of mechanical non-integrable system are broken, thus a chaotic 
motion in phase space takes place. Chaos, most likely with mixing properties, 
is very frequent in mechanical systems. [21], [2 2] 

A very important and popular example of uniformly mixing transforma- 
tion is the, so called, baker transformation that operates in the phase space 
X = [0, 1] X [0, 1] and it is defined by the following procedure: 

i. -squeeze the 1x1 square to a 2 x | rectangle, and, 

ii. -cut the rectangle vertically into 2 rectangles and pile them up to form 
another 1x1 rectangle. 

in doing so the point of the square will move as: 

{x,y)^S{x,y) = l ^^^''f^^'i ^ '{■■'i^'^^f (4.3.2) 

I (2a; - 1, i |y), if : ..^ < x < 1 ^ ' 

The transformation is shown in fig. 1, where in the first square, the one 
corresponding to i = 0, the lower half is shadowed and corresponds to a 
subset B. It is easy to see that this transformation is reversible. The fate of 
this area B, evolving to the future, is shown in the right side of the figure, (it 
is transformed in a great number of horizontal strips with area 1/2) and 
evolving to the past, in the left side (the strips are now vertical). An smaller 
subset A is also shown. It is then easy to verified that condition (4.3.1.) 
is fulfilled (the final measure of StiB) fl A will be \yi,i,{Ai)^ being the initial 
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measure of B just | so eq.(4.3.1) is satisfy). We will study this transformation 
in all detail in the subsection 4.5 using the fine-graining technic. 

Much more complicated mixing evolutions than (4.3.2) can be invented. 
In fact, baker's transformation is the simplest of all, it is the simplest model of 
the famous Gibbs ink drop. Gibbs try to explain the essence of irreversibility 
with the ink drop model. If a drop of blue ink is introduced in a glass of 
water, even if the volume of the ink drop remains constant (as the volume 
of any subset of mechanical phase space, according to Liouville theorem) we 
will have, after a while, an homogeneous mixture of bluish water. What it 
is happens is that the motion of the water is mixing and therefore the ink 
drop is deformed (even if its volume is constant) in such a way that it is 
transformed in a set of very thin filaments that are present in every portion 
of the water giving the sensation that the water has become bluish. The 
growing of this filaments-like structure gives an arrow of time and it is for 
Gibbs the essence of irreversibility. This phenomenon is modeled by the 
baker's transformation. In fact, let us consider a small rectangle axb within 
the square 1 x 1, let say a small task of lower quality fiower within the bread 
mass. The height of the task will successively became: jb, ...while 

the base of the task will became: 2a, 4a, ...ta, ....in such a way that the area 
is conserved. Eventually a time will arrive such that ta > 1 and then the 
task will be cutted in two, and then in four, eight, etc., in such a way that 
it will become a set of horizontal filaments of decreasing height, namely a 
"cubistic" picture of the ink drop, so bakers transformation is just a model 
of the ink drop phenomenon. 

If now we consider the much more complicated evolution of the ink drop, 
and if the volume of the ink drop is the 1% of the volume of the water, it is 
clear that the motion of usual water is mixing according to definition (4.3.1), 
as bakers transformation. In fact, if the motion is mixing, when i — > oo every 
subset A C X will have a 1% of ink and, therefore, the distribution of ink 
will become homogeneous. As this is the case with the real ink drop we can 
conclude that the real motion is mixing. 

It is a straightforward consequence of the definition that p*-mixing implies 
p*-ergodicity. In fact, if B is an invariant set eq. (4.3.1) reads: 

li^{AnB) = ^x^{A)^x^{B) 

for all set A. Now if we take B — A we obtain ijl*{B) — [//*(5)]^ and 
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therefore either ii*{B) = 0, or ii*{B) = 1, so ii^{X — B) = 0. So the evolution 
is p*-ergodic. 

Now we have arrived to our most important 

Theorem 4.3.1. [6].- Let St be an ergodic transformation, with stationary 
density of the associated Frobenius- Perron operator, operating in a phase 
space of finite p*-measure. Then St is mixing if and only if {Ptp} is weakly 
convergent to p*, i.e., 

lim {Ptp,a)^{p,,a) (4.3.3) 

t— »oo 

where cr is a bounded measurable function. 

If a sequence is weakly convergent it is also Cesaro convergent, so we can 
see again that mixing evolutions are ergodic. 

So the mixing property assures a weak convergence of {Ptp} to p*. But, 
e. g., in the example of the baker transformation the strong limits toward the 
far past or the far future do not exit. In fact, the support of any distribution 
function (if it has a measure < 1) will be a set of infinity horizontal toward 
the future or vertical lines toward the past. These set can not be the support 
of any regular distribution function. Nevertheless the weak hmit (4.3.3) do 
exists with p* = 1. 

The physical meaning of theorem 4.3.1 is very clear: Let us consider a 
(non-viscous) fluid in motion in a cubic box. As energy is conserved the 
motion will never stops, and therefore, according to the laws of mechanics, 
equilibrium will never be attain, and Ptp will have no limit. This will be the 
case if the motion is oscillatory, namely a pressure wave that oscillate back 
and forth between two parallel walls of the box. But if the motion is mixing it 
is so complicated that there are portions of the fluid moving in every direction 
near every point of the box. In this case if we take the inner product {Ptp\a) 
we are making an average that goes to an equilibrium average (p*|cr) when 
t — >■ oo. Therefore, even if there is always motion, the motion average gives 
an image of equilibrium. This is the profound meaning of theorem 4.3.1. and 
the way to obtain a synthesis of the apparent contradiction of dynamics and 
thermodynamics: 

-even if the dynamics says that the energy is conserved and the motion 
will never stops 

-there is a thermodynamical equilibrium in average, because the motion 
is mixing. 
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From this point on fine-graining and coarse graining follow different paths, 
as we have explained in the introduction and we will discuss below. But let 
us remember that the problem is not completely solved, since all the nice 
inequalities of subsection 3.6, that are necessary to explain the second law , 
are equalities for reversible system and all system in nature are considered 
to be, al least microscopically, reversible. 

There are system endowed with properties more chaotic than mixing, 
they are: 

i. - Kolmogorov systems [6] that necessarily are mixing [23]. 

ii . - Anosov systems- [6] , [22] , [24] . 

iii. -Bernoulli systems, the most chaotics of all [25] . Baker transformation 
is, in fact, a Bernoulli system [26].. 

4.4 Exactness. 

We will now introduce a property that (apparently) will solve all our problems 
If St is a p*-measure preserving transformation operating in a phase space 
X, then St is said to be p*-exact if: 

lim i2,{St{A)) = 1 (4.4.1) 

I— >oo 

for all sets A of non zero measure. This is possible even if St is p*-measure 
preserving since a evolution is measure preserving if eq. (4.1.2) is satisfied 
and this equation is not equivalent to iJ,^{St{A)) = /x*(A), if the evolution is 
not reversible. Renyi map is a good example. 
Let us consider a dyadic Renyi map: 

R : [0, 1) ^ [0, 1), ...X ^Rx = 2x{modl) (4.4.2) 

As the length of any subset A is multiplied by two in each transformation, 
this map is exact since it satisfies eq. (4.4.1). Anyhow it is also measure 
preserving. In fact, let as consider, e. g. the subset A — [0, |), R~^{A) is 
[0, |) U [|. |) and. therefore, both subsets have measure |. 

If p* is the uniform density of the microcanonical ensemble we say that 
St is uniformly exact. 

The essential think to understand is that reversible system cannot be 
exact. In fact, for reversible p*-nieasure preserving transformation we have: 

MSt{A)) = i^4Si\St{A))] = ii^A) (4.4.3) 
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thus the definition of exactness is violated. Since usually classical dynamical 
system are measure preserving, by the Liouville theorem, and reversible they 
are not exact. Nevertheless, as we shall see exactness is really the property 
we are looking for. Precisely: 

Theorem 4.4.1. [6] If St is a p^,-measure preserving transformation op- 
erating on a finite normalizablc phase space X and Pf is the associated 
Frobenius- Perron operator corresponding to St, then St is p*-exact if and 
only if: 

lim|| Ptp-p* 11=0 (4.4.3) 

t— *oo 

Therefore: Ergodicity corresponds to Cesaro convergence. Mixing corre- 
sponds to weak convergence, and, exactness corresponds to strong conver- 
gence (i.e. convergence in the norm). A strongly convergence sequence is also 
weakly convergent, thus we can deduce that exact evolution era also mixing 
evolution and therefore ergodic evolutions. Moreover, since we are looking 
for a strong limit we see that working with ordinary distribution functions we 
will find this limit only if the transformation is exact, but ordinary classical 
(microscopical) system arc not exact since they are reversible and measure 
preserving. As an example we have shown that the reversible baker transfor- 
mation has not strong limits toward the past and the future, in fact baker 
transformation, being reversible cannot be exact. Thus our problem is now 
clearly stated: if we want a strong limit our evolutions must be exact, but 
exact evolutions are not reversible and all microscopical transformation are 
reversible, therefore we can not have a strong limit. Furthermore we have 
also the 

Theorem 4.4.2. [6] .Let Pt be a Markov operator operating in phase space 
X. Then the conditional entropy of Ptp with respect to density p* goes to a 
maximum value of zero as i — > oo, i.e., 

lim Hc{Ptp\p.) = (4.4.5) 

t — >oo 

if and only if Pt is p*-exact. 

This theorem tell us the necessary and sufficient criteria to be able to 
state the second law of thermodynamics in third-order form, namely, for the 
entropy of the system to converge to its maximum value regardless of the 
way in which the system was prepared. This condition is that the system 
must evolve according to an exact transformation. But such systems do 
not exist in nature. So dynamics cannot be related, at least trivially, with 
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thermodynamics. Therefore our theory must be modified one way or the 
other. 



4.5 Mixing studied by the fine-graining technic. 

Mixing evolution are studied by the fine-graining technic in papers [26], [27], 
and [28], using a perturbative method, that can be implemented in any ex- 
ample. For didactical reason we will present the most important mixing 
evolutions and we refer to the papers above for the general perturbative 
method. 

4.5.1 The Renyi maps. 

The /5— adic Rcnyi map R on the interval [0, 1) is the multiplication, modulo 
1, by the integer /3 > 2 : 



The forward iteration of the Renyi map n times, define a " cascade" or time 
evolution with time t — n & Z. This evolution preserve only the Lebesgue 
measure, as we have shown after eq. (4.4.3). The density functions p{x) 
evolve according the Probenius-Perron operator U : 



Gcl'fand-Maurin theorem 4.A.I. tell us that we can found an spectral ex- 
pansion in the eigenvectors of this operator in an adequate rigged Hilbert 
space. In fact, using the perturbative methods of papers [26], [27], and [28], 
the spectral decomposition of U can be found and reads: 



R : [0.1) [0, 1) : ....x ^ Rx ^ Px..{mod..l), 



(4.5.1.1) 




(4.5.1.2) 



OO -1 ^ 

n=0 t-' 



•n 



(4.5.1.3) 



where : 



n-l 



S„(x)) = |x"+ ^ 



X' 



.m 



m!(n — m)! 




(4.5.1.4) 



m=0 
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where the Bn{x) is the n-degree Bernoulh polynomial defined by the gener- 
ating function: 

^=£m» (4.5.1.4) 

- 1 ±1 n! ^ ' 



and: 



~ (1| n = 



n! 



where (1| is the constant distribution function. From eq. (4.5.1.5) we can see 
that the elements of the spectral decomposition (4.5.1.3) do not belong to C 
but to a larger space where the Dirac 5 must have a precise mathematical 
meaning. This space is, in fact, a rigged Hilbert space that we shall define 
below, in full agreement with the Gel'fand-Maurin theorem. The system 
{\Bn)i {Bn\} is bi-orthonormal and complete, namely: 

(S„|SJ = 5„^ (4.5.1.6) 

oo 

Y,\Bn){Bn\^l (4.5.1.7) 



n=0 



The spectral decomposition (4.5.1.3) acquires a precise mathematical mean- 
ing if we define, as test space $ the space of polynomials V .This space is 
dense in L—L? , nuclear (in fact, it is the union of an infinite and discrete 
set of finite dimensional spaces), complete, stable under U . and U is contin- 
uous in the topology of V . It is, therefore an appropriate test space to give 
a meaning to the spectral decomposition whose elements belong to But 
other kind of test functions spaces can be defined and we will obtain differ- 
ent spectra, e. g.: a continuous set of eigenfunctions can be found, with and 
adequate rigging, showing that the Renyi map have continuous spectrum, 
precisely the set of complex numbers z such that |2;| < 1 (all the mixing op- 
erator have an spectral decomposition with a continuous spectrum in Hilbert 
space). 

liteZ from eqs. (4.5.1.3,6) we can see that the evolution operator is: 

oo 1 ^ oo -1 _ 

U'-T. ^t\Bn){Bn\ = |1)(1| + E -j|5„)(S„|....(4.5.1.8) 

n=0 n=l 
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If we would like to work in space L only we must remember that all the 
formulae above are just weak equation, e. g. the last one is just: 

CXD 1 ^ 

(<7|C/V) = (^|1)(1|P) + (^1 E ^\^n){Bn\p) (4.5.1.9) 

n=l 

for all p,(7 e and p is a density. Then as /3 > 1 and p is normalized we 
have: 

lim (cT,[/V) = (cx|l) (4.5.1.10) 

in perfect agreement with theorem 4.3.1. Eq. (4.5.1.9) is just the weak 
version (or coarse graining) of eq. (4.5.1.8) that allows us to work within 
space C but using always weak limits as (4.5.1.10). But if we work with 
functional directly, namely in space from eq. (4.5.1.8) we can say that: 

lim C/*|p) = 11) (4.5.1.11) 

which is a strong limit, namely the fine graining version of (4.5.1.10). If we 
call /3 = < 1 , from eq. (4.5.1.9) we can also say: 

p{t) = U'p = + pi{t)e-^', ...p* = |1) (4.5.1.12) 

where = |1) is the equilibrium distribution function and pi{t)e~'^ is some- 
thing like a "fluctuation" around the equihbrium state. We write this last 
equation because we will find a similar equation in the quantum case. 

4.5.2 The baker's transformation. 

The /3-adic = 2, 3, .. baker's transformation in the unit square Y — [0, 1) x 
[0, 1) is a two-step operation: 

i. -squeeze the 1x1 square to a /? x ^ rectangle, and, 

ii. -cut the rectangle vertically into p rectangles and pile them up to form 
another 1x1 square. 

Then: 

(x, y) B{x, y) = {f3x-r, ^4—)' --for ■^<x< ..r = 0, .../3-1 (4.5.2.1) 



47 



This equation is an obvious generalization of eq. (4.3.2), which is the par- 
ticular case of eq. (4.5.2.1) for /5 = 2. As we can see we have a sort of two 
Renyi maps one in each coordinate. Baker's transformation is a Bernoulli 
shift and has Kolmogorov Sinai entropy log 2/5 [25]. The invariant measure 
is Lebesgue measure. The density function p{x, y) evolves according to the 
Frobenius-Perron operator U : 

Up{x, y) = p{B-\x, y)) = p{^, f3y-r), ..far ■^<y< -r = 0, ..(4.5.2.2) 

This operator is unitary in the Hilbert space £ = the equilibrium distri- 
bution function is the constant function = 1, and the Lebesgue spectrum 
is the unit circle plus the simple eigenvalue 1 

As the baker's transformation B is the natural extension of the Renyi 
map R the conclusion that we can obtain are the same and we refer to [29] 
for details. B acts on the Liouville-Hilbert space C — L"^ = L^. ® L"^ and 
a suitable initial biorthonormal system can be constructed from the tensor 
products of the eigenfunctions of the /3-adic Renyi map (cf. eq. ((4.5.1.4) 
and (4.5.1.5)) 

Wnm >= Br,{x) Bm (y)-and... <<Prim \ =-Bn {x)Bm{y) (4.5.2.3) 

Using these bases and the same perturbative method as before the following 
spectral decomposition can be obtained: 

U = l/oo X/oo l + E E 1^\U ><f.,r I + E ><U\ ...(4.5.2.4) 

u=l U=0 r=0 ) 

where the vectors \f^^r > and <f^r I ^e obtained from the vectors of 
eqs. (4.5.2.3). As we have said the Liouville spectrum is the unit circle plus 
the eigenvalue 1, so in the new spectral decomposition we have found new 
eigenvalues 1/(3" < 1. 

The initial vectors (pnm and are linear functionals over the spaces 
= (g) and = Vx 'Si Ly. Furthermore it can be shown that the 
vectors fnm £ and fnm e ^+ are also functional over the same spaces, 
so the spectral decomposition (4.5.2.4) can be implemented if we use these 
functional vector spaces. We have enlarged the state space with densities 
that can be distributions in the y coordinate, in the case of e. g.: if the 
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y distributions are Dirac's deltas we will have a distribution whose support 
is a set of horizontal straight lines, that we shall call a "horizontal Dirac's 
comb" . In the case of we must change the y by the x and we would have, 
e. g.:" vertical Dirac's combs" 

Now, "mutatis mutandi", we can repeat what we have said in eqs. (4.5.1.9) 
to (4.5.1.12), and we will find similar equations for the baker's transforma- 
tion. The equilibrium distribution, in this case will be: 

P* = l/oo > (4.5.2.4) 

where |/oo >e 

4.6 Appendix 4. A. Rigged Hilbert spaces [10], [29]. 

As it is well known all linear spaces of the same dimension are isomorphic 
if this dimension is finite. This is not the case if the dimension is infinite. 
In fact, let us consider the infinite sequence of the vectors of a basis of an 
infinite dimension vector space 

{$„:n=l,2,...}, (4.A1) 

Let V be the vector space of all finite linear combination of the vectors of 
the basis above, namely: ^' e V if: 

* = X:cA (4.A2) 

n=l 

y is a linear space of infinite dimension, but we will see that we can built 
other spaces using basis For instance we can add to V the limit points 

of all the convergent infinite sequences of vectors of V. But defining different 

criteria of convergence wc will have different set of limit points and, therefore, 
different vector spaces. The most use full convergence is the convergence in 
the norm. The sequence {^i} converge in the norm to a limit point % if: 

lim|| X - *i Ih (4.A.3) 

j— >oo 

If the sequences {^i} are sequences of vectors of V and they converge in 
the norm and we add the limit points x of these sequences to V we obtain 
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a larger space Ti where we have finite sequences like (4. A. 2) but also limit 
points of infinite sequences. We will say that Ti is the closure of V and also 
that V is dense in H.. But we can use other kind of convergencies, namely 
other topologies, and we will obtain different spaces. 

Let us suppose that we chose the sequences, such that the coefficients c„ 
satisfy the condition: 

oo 

E Icnr < OO (4.A4) 

n=l 

Then, adding the corresponding limit points, we obtain a Hilbert space Ti. 
which contains all the sequences that converge in norm, and it is call also 
the competitions of V respect to the topology of the norm. But we can also 
consider a the infinite dimension linear space S of all, either finite or infinite, 
linear combinations of the basis {$n}; namely all the linear combinations 
^ = Cn^n with limitations impose over the coefficients c„. Of course we 
cannot define a norm in such a space, but we now we have three infinite 
dimensional linear spaces such that: 

y C C S (4.^.5) 

Let us define the inner product: 

{Lh)-Y.Kcn (4.A6) 

n 

Then Ti is the space of the vectors h = X^n Cn$n such that {h,h) = I] |c„p < 
oo. Let us now define the conjugated space of 7i, Tt^ C S of all linear 
functional over Ti namely vectors / = J2n bn^n such that the inner product: 

/W = (/,^) = E^:cn (4.A7) 

n 

is convergent for all h E H. The convergence of this inner product is a 
consequence of Schwarz inequality: 

\{f,h)\'<{h,h){fj) (4.A8) 

so (4. A. 7) converge if (/, /) = l^nP converge and, therefore, / G so 
Ti = Ti^. (We can as well define the space of antilinear as f[h] = {h, f ) i.e. 
bra are linear functional and ket can be consider like antilinear functional) 
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Let us now define a new space ^2 as the space of all vectors uj = J2n "^n'&ra, 
endowed with coefficients c„ such that-they satisfy the following set of infinite 
conditions: 

J2 lunl'^rT < oo, ...m = 1, 2, 3 (4.^.9) 

n 

Obviously Vt (ZTi. Let us now find the conjugate space of fi, Vt^ C S, namely 
the space of convergent linear continuous functional over VL. These functional 
read a — Y.n'^n^n and they are such that: 

a[uj\ = (a,cu) = ^<M„, (4. A 10) 

n 

is convergent for all a; £ O. Therefore: 

Mn-"^ < oo, ...m = 1, 2, 3...(4.A11) 

n 

In fact: according to Schwarz lemma we have: 

|^<n-T«„nT|2< iY,\v^\'n-A lu^^n-) < +oo, (4.A12) 

and the r.h.s. is convergent if eqs (4. A. 9) and (4. A. 11) are fulfilled. 

As it is obvious that — S, we now have the following set of infinite 
dimensional spaces: 

V ^^1^71 = ^ cVf" C = S (4.A13) 

Any triplet: 

QgHg (4.^.14) 

like those of eq. (4. A. 12) and others that can be obtained, e. g.: if limit 

the m's of eq. (4. A. 9) to be just 1< m < M, for some M & N, arc called 
Gcl'fand triplets. Q is known as the test space and , as the rigged space. 
Mathematically it is convenient that the test space would be a nuclear space. 
Heuristically speaking nuclear spaces are the infinite dimension spaces en- 
dowed with the largest number of properties of finite dimensional spaces, 
among then they have discrete spectral decomposition. Precisely, nuclear 
spaces are spaces obtained, so to say, as the union of an infinite sequence of 
spaces of finite dimension. As space fl, from our point of view, is the space 
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of operators corresponding to real measurement apparatuses and as these 
devices make only a finite (so less than discrete) number of measurements; 
logically must be a nuclear space. In fact., even if physical devices make 
a finite number of measurements, we can conceive that this numbers grows 
with the progress of technology. Then an infinite, but discrete, number of 
measurements would corresponds to the limit of an infinitely long period of 
technological progress. A finite number of measurements will corresponds to 
a test space of a finite number of dimensions. Then the test space correspond- 
ing to the limit of technological progress will be a nuclear space, since this 
space is the limit of a sequence of finite dimensional spaces. E. g.: a measure- 
ment device make n measurements, that can define n points of a curved, that 
can be interpolated by a polynomial of degree n. The space of polynomials of 
degree n will be the test space that corresponds to this device. In the limit 
of technological progress the test space will be the space V of polonomial 
of any degree, in fact a nuclear space. Generally speaking choosing different 
nuclear test functions spaces we can also chose the physical properties of our 
measurement devices. 

In finite dimensional vector spaces the eigenvalue problem, for every 
selfadjoint linear operator A can be solve in a unique way. Namely we 
can find a unique spectrum {flj} and an orthonormal basis {^n} such that 
A^n = fln^n-This is uot SO for infinite dimensional linear spaces, since the 
spectrum depends on the rigging we use, nevertheless it can be demonstrated 
the Gel'fand-Maurin 

Theorem 4.A.I. If ^4 is a self adjoint operator in TC there is always a 
complete set of eigenvectors of A in some rigged Hilbert space fl^ . 

Let us give to very important examples: 

i.-Let S be the space of functions f{x) of one real variable x and let 
A — P — —i-^ he the self adjoint momentum operator in H — L^. The 
eigenvectors of P are the plane waves e^'^^ with do not belong to since 
they have not finite norm. Nevertheless they can be considered as functionals 
over a convenient space test function $ since: 

/+0O 
e-'''f{x)dx^f{k), (4.A15) 
-oo 

where f{k) is the Fourier-transform of f{x) and $ is any subspace of Ti. 
such that cq. (4.A.15) id convergent. Then we have the Gel'fand triplet 

^ CH C^"^ and e*''""' G 
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ii.-Let S be as in the example above and A = Q = x he the position 
selfadjoint operator in 7i = L^. The eigenvectors of Q arc the Dirac's deltas 
5{x — y), since QS{x — y) — yS{x — y) , these distributions do not belong to 

since they are not even functions. Nevertheless they can be considered as 
functionals over a convenient space of test functions since we can rigorously 
define these deltas as the functionals: 



where f{x) is any function of $. Usually physicist write this last equation as: 



Syifix)] = {Six - y), fix)) = Six- y)fix)dx = fiy), (4.A.17) 



even if the integral in this last equation have not a rigorous definition. Usu- 
ally $ is the set of function with nice properties e.g. they are continuous, 
derivable, with compact support, etc. Then dy G 

these examples show that usual operator don not have their eigenvalues 
in Ti but in properly chosen rigged Hilbert spaces. 

5 The Quantum Evolution. 

As the laws of quantum evolution are well known (cf. [10], [17], [18]), in this 
section we will see the use of the no-graining and coarse-graining technics in 
quantum mechanics. 

5.1 The case of discrete spectrum. 

Let us begin making just an heuristic calculation. Let Ti. be the quantum 
Hilbert space. Let {\i >} be a energy eigen-basis of this Hilbert space, where 
i is a discrete index. The quantum Liouville space is C — Tl x H, and a 
generic density matrix reads: 



(4.^.16) 






where since p = p\ it ispij = p*^. 
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Let O be a self adjoin operator, it reads: 

O^J^Oij\i><j\ (5.1.2) 

where Oij — O*^. 

The mean value of operator O.in the quantum state p is: 

< O >p= {p\0) = tr{p\0) = J2PijOji (5.1.3) 

ij 

as |i > is an energy eigen state we have: 

H\i >= uji\i > (5.1.4) 

where uji is the energy of state |i >. The time evolution of this eigen states 
reads: 

\i(t) >= e-*^'*|i > (5.1.5) 

Therefore the time evolution of p is: 

p{t) = EijPiMt) >< jit)\ = Ei,- A,e^(-^-"^-)*N X J\ 
= EiPiili >< i| +Ei^,-Pi,-e*^"^-"^)*|i X (5.1.6) 

Then the time evolution of the mean value of eq.(5.1.3) is: 

< O >p(t)= (p(t)|0) = YlpiiOii + E (5.1.7) 

Now let us suppose that the steps of the spectrum arc so small and the 
function under the second summatory of the r.h.s. of the last equation is so 
nice that this summatory can be approximated by an integral. Therefore, if 
the function is nice enough, from Riemann-Lebesgue theorem we would have: 

lim {p{t)\0) = Y.P^^O^^ = (p*|0) (5.1.8) 

t— »oo ■^"^ 

I 

where we have defined an equilibrium density matrix p^ij = puSij. This 
equation would be the quantum analog of the classical equation (4.3.3) for 
mixing system and it would show that both system have a similar behavior 
and opens the possibility to use classical theorems in the quantum case also. 
Of course this demonstration is not rigorous, but it serve to motivate the 
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study of continuous spectra of next section. Using continuous spectra we 
will find a rigorous theorem. The role played by continuous spectra in this 
case is not strange since evolution operators of mixing systems have this 
kind of spectrum. [30]. Anyhow, we can also say that the most we can get 
is a weak limit, since Riemann-Lebesgue theorem cannot be used directly in 
eq.(5.1.6). Furthermore, de decomposition of the r.h.s. of this equation is 
not a decomposition within space Ti since its second term has null trace (cf. 
eq.(2.3.2)). 

5.2 The case of a continuous spectrum. [31] 

In the next subsection we will consider the Priedrichs model, with can be 
defined in Hilbert space H with a energy eigen-basis {|1 >,\oj >},0 < u < 
oo, with hamiltonian operator: 

H = uji\l >< 1| + /o°° dujuj\(jj X uj\ 
+\J^ dujg{u)[\u >< 1| + |1 ><uj\] (5.2.1) 

In this section this formula will be only used as an example of an operator 
expanded in a continuous spectrum basis to conclude that the expansion of 
a generic self adjoint operator reads: 

roo f- /'OO 

0=1 dwO^luj >< uj\+ dujduj'0^^>\u> >< uj'\ (5.2.2) 

Jo J Jo 

where O^;, Oj^w'^re regular functions such that O^j G R,0*,^^ = O^joj'- Below 
we will say that O G $, a space with some properties that we will chose for 
convenience. Thus functions Oi^O^juj' will be restricted by this choice. 

The first term of the r.h.s. of eq. (5.2.2) will be called the singular 
component of O , since it could be written as the second term but with a 
singular coefficient Ou^oj' — OujS{u; — a;'). The second term will be called the 
regular term. 

Let us consider density matrix at time t = 0: 

P POO 

p(0) = J J duduj'p^^'luj Xuj'l (5.2.3) 

at time t this state reads: 

p{t) = J dujduj'p^^>\uj >< uj'\e-'^'^-'^'^^ (5.2.4) 
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If we consider that O can be written as the p of eq. (5.2.3) but with 
coefficients Ouj5{uj — ui') + 0^^^^l (as in eq. (5.2.2)) the mean value of operator 
O in the state p{t) is: 

< O >,(i)= {p{t)\0) = trip{t)0) = !^ dup^^O^ 
+ //o"rfa;rfa;'p,,,0,,,e-^(--')*.... (5.2.5) 

Now,if space $ is chosen in such a way that Riemann-Lebesgue theorem can 
be used, namely let the functions of $ be L^, we have: 

roo 

\\-m<0 >pit)= I dujO^p^^ (5.2.6) 

As this equation is valid for any operator O e $ we may try to find a density 
matrix such that: 

lim< O >p(t)=lim {p{t)\0) = (p,|0) (5.2.7) 

t— *0O ^ I— >oo 

It is easy to see that the density matrix cannot be find if p^^/ is a regular 
function of variables uj,uj', i.e. from eq. (5.2.5) we see that to obtain this 
result it is necessary that p^^t^/ = 0, uj ^ uj' and p^^^^ ^ 0, but we cannot write 
Pww' = Puj^ij^—'^')-! because in this case the p is not regular. [32]. Then we are 
forced to consider states with diagonal singularities, namely with the same 
operator's pathology. So we are forced to introduce singular components in 
the density matrix, namely thepaj,a;' of eq (5.2.3) cannot be regular and it 
must read something like Puj^iuj — uj') + Puj,uj'- But now if we try to find the 
mean value (5.2.5), the Ouj5{u) — cu') term and the PojS{u! — uj') term produce 
the result: 



p poo poo 

/ / OuiS{uj — uj')pi^i5{uj — ijj')d(jjd(jj' = / Oi^pi^6{Q)d 
J Jo Jo 



;a; — >• oo 



which is divergent. Therefore to have a formalism free of these problems we 
are forced to make a fresh start and to consider that the operators O are 
defined by the regular functions Oi^,0^^^' and the state functions p are the 
matrices of rigged space $^ defined as the linear operators on space $ and 
therefore are defined by two regular functions Poj,Poj,oj'- Then we have: 

poo p poo 

(p|0) = dwpu,0^ + J duduj' p^oj'O^uj' (5.2.8) 
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where G R,plji^ = Puu'- So we are forced to introduce a singular compo- 
nent Pi^ also in the density matrices. Now can be found, it is the functional 
of space with ^ O, p^^^i — 0. The consistence of this method is proved 
by the logical physical results of paper [31]. 

Eq. (5.2.7) can now be consider as the rigorous quantum analog of the 
classical cq. (4.3.3). Wc can call the weak limit of this equation the "quantum 
mixing" property and state the following 

Theorem 5. 2.1. -Quantum system with continuous spectrum are en- 
dowed with the quantum mixing property, (provided we use the formalism 
based in eq. (5.2.8)) 

Eq. (5.2.7) can also be considered a prove of a weak dccoherence in 
quantum systems. This would the no-graining conclusion. But we would like 
to have a strong decoherence. Then we can follow two ways: We can use 
coarse-graining. This technic is well known, so we refer to papers [33] and 
[34]. Or we can use fine-graining. In this case, in order to obtain a strong 
limit from the weak limit of eq. (5.2.8) wc must gave a precise sense to all 
terms of the r.h.s.of eq. (5.2.5), rigging the Hilbert space Ti in such a way 
that all the mathematical characters are well defined. So working with the 
functional of space we can write the strong limit: 

lim p{t) =p, (5.2.9) 

t— »oo 

that corresponds to the classical strong limit (4.4.3). 

Now we would like to obtain, not just a limit, but the time irreversible 
evolution of p{t) that yields the limit (5.8.9). Unfortunately our silkiness to 
work with continuous spectra is very limited [35] so we are forced to use some 
mixed technics, as we shall see in the next subsections. 

5.3 Friedrichs model. 
5.3.1 The general formaUsm. 

We believe that the, well known, Friedrichs model is the best quantum exam- 
ple to fix the ideas. In this example we have a free (naked) stable quantum 
state |1 > (which is postulate to be real: K\l >— |1 >) that becomes un- 
stable when coupled to a continuous field \uj > . The stable state can be 
considered as a simplified model of an atom in an exited stable state, that 
becomes unstable if coupled with an electromagnetic field, which, in the 
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model, is represented by the continuous field. Thus, let us consider a Hilber 
space 7i, with a basis {|1 >,\uj >}, < u < oo such that: 

< 1|1 >= 1, < >= 0, < u;\u;' >= S{u; - uj')....{5.3.1) 

roo 

1 = |1><1|+/ \uxuj\du (5.3.2) 

and a system with free hamiltonian: 

/■oo 

Ho = uji\l><l\+ a; I a; >< a; Ilia; (5.3.3) 

Jo 

and cui > 0. Therefore the spectrum of Hq is with a degeneration at cui. 
Let the interaction hamiltonian be: 

roo 

Hi = X g(to)(\l >< u}\ + la; >< l\)du; (5.3.4) 

Jo 

where g{u!) is an interaction function endowed with all sort of nice properties: 
it is analytical, well behaved at a; — > +oo, etc., etc. The total hamiltonian 
is: 

H = Ho + Hi (5.3.5) 

This hamiltonian can be diagonalized, using standard technics. Then we 
obtain: 

adv adv 

H= u!\u!,ret>< LJ,ret \du! (5.3.5) 

Jo 

adv 

where {\u},ret>} are the usual retarded or advanced bases [38]. We can see, 
comparing eq. (5.3.3) and eq. (5.36) that the interaction has erased the 
discrete component of the spectrum. In fact, state |1 >has became unstable 
and now it is just a pole in the corresponding S-matrix. Any how using 
eq. (5.3.6) we can compute the time evolution of any state, e. g. the state 
|1 > at i 0. As we have just said state |1 > of the free system (5.3.2) is 
transformed in an unstable state by the interaction (5.3.4), in such a way 
that the survival probability P(t) = | < l|l(t) > p vanishes when t +oo. 
It is also known that P{t) has a vanishing derivative when t = (Zeno effect), 
then it has a decreasing exponential behavior, and finally oscillates for big t 
(Khalfin effect) (fig. 2) [36]. 
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5.3.2 Hilbert and rigged spaces. 

Let us forget for a moment the problem of the unification of the dynamics 
and the thermodynamics and let introduce some equations related with the 
problem of the time asymmetry, as it was stated is subsection 1.1. As we 
shall see the previous equations are all what we need to define the quantum 
arrow of time according to the coarse graining school.. As this school always 
work within Hilbert space Ti. , the following property hold: 

K-.n^n (5.3.7) 

Instead for the fine greaning school we need a richer structure. In fact, we 
need to define two subspaces 0± C 7i , in a convenient way. To do so let us 
consider a vector >e 7i and its components < uj\ip > and let us promote 
the real energy a; to a complex variable then: 

Iv? >G 0± iif < z\^>e Hlf\S (5.3.8) 

where are the Hardy classes from above and below respectively (cf. App. 
5. A) and S the Schwarz class of functions. It can be roved that 0± are nuclear 
spaces. Then we can then define two Gel'fand triplets: 

C H C0!! (5.3.9) 

(p+GHG (p'l (5.3.10) 

5.3.3 The rigged Hilbert space formalism. 

Using analytical continuation technics (cf. [15], [36], [37], [38]), essentially just 
the Cauchy theorem, it is possible to obtain a new spectral decomposition of 
the identity operator 1 and the hamiltonian operator H as: 

1 = \zi, - >< ^1, +1 + J \z, - >< z, +\dz (5.3.11) 

H = zi\zi, - >< zi, +\+ j - X (5.3.12) 

where \zi, — (p^, \zi,+ (p^ , is a complex rot of equation a{z) — 
where: 

a{z) ^z-u;, + X [ ^*^^*^^^^^ dz (5.3.13) 

Jr z — uji 
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and r is any curve that goes from the origin of the complex plane to the 
positive infinity of the real axis and passes below zi (fig. 3). The first 
terms of eqs. (5.3.11) and (5.3.12) are produced by the residues of the poles 
corresponding to the roots located at the zeros of equation a{z) = or, what 
is the same thing, the poles of the S-matrix. We can see that the discrete 
component of the spectrum, that we have lost in eq. (5.3.6), reappears in eq 
(5.3.12) in the form of a matrix of the rigged Hilbert space. 

Now we have several possibilities to chose the curve F, that are used by 
different authors: 

i. -To use a generic curve F. 

ii. -To use curve F' of fig. 4, in such a way that, as the vertical paths of 
the curve are mutually cancel we are mostly integrating on the real positive 
axis. 

iii. -To take the negative real axis as the integration path. 

iv. -To define a tilde operation as: 

f{x)g{x)dx 

for all g{x) in the test function space. In this case the complex integral 
formally becomes a real one. 

If we use this last method and forget the tilde eq. (5.3.12) reads: 

/■oo 

H = zi\zi, - >< zi, +1 + / u!\u!, - >< (jJ, +\duj (5.3.12') 

Jo 

so we have build a basis {|2;i, — >, — >} for space 0_(cf. [37] for details). 
These vectors reads: 

\zi, - >=< llzi, - > (11 > +A / dou ^ \ \uj >) (5.3.12") 

Jo [Zi — U!\- 

\u,->^ \u > +M^(|1 > +A rdu'^^lu' »....(5.3.12"') 

where the subindex "-" in the denominator of the integral in eq. (5.3.12") 
means that the curve F' must be used for the integration. 

Now we have two spectra to compare: (5.3.5) and (5.3.12"). The main 
difference is that (5.3.5) is structurally unstable when A ^ 0, while (5.3.12") 
is stable. In fact, an algorithm is call structurally stable if it does not change 
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much under small changes of the coefficients. When A = the spectral 
decomposition of H is (5.3.3). If A is small, a small change of A that makes 
A = 0, produces a big change in the usual decomposition, that goes from 
(5.3.5), with no discrete term, to (5.3.3) with the discrete term uji\l >< 1| . 
The sudden vanishing of this term when A — > is a catastrophe (precisely a 
Poincarc catastrophe) that creates a great number of problems if we try to 
perform an expansion around A = in Hilbert (or Liouville) space. On the 
contrary (5.3.12") is stable, since it has the term zi\zi, — >< zi, +|that goes 
to a;i|l >< 1| when A — as we shall see. 

From eq. (5.3.12) it can be seen that — > and < -Zi, +| are respectively 
the left and right- eigenvector of H corresponding both to the eigenvalue zi_ 
It can be proved that: 

< >= 1 

< z + \z,- >= 

< Z, +1-21 _ >= 

<z. + \z, -' >= d{z - z') (5.3.14) 

It can also be proved that: 

< zi, -\zi_ >= 

< zi, + \zi, + >= (5.3.14') 

namely there are non-null vectors of zero norm in spaces 0^ and 0^ [39]. 

Let us call Zi = Pi — |7i, where 71 > 0. Then from eq.(5.3.12') we can 
obtain the time evolution of \zi{t),— > and \zi(t),+ > precisely: 

\zi{t), - >= e-''''\zi{0), - >= e-^^i*e-^*|zi(0), - > 

\zi{t), + >= e-'"i*|zi(0), + >= e-"^''e^'\zi{0), + > (5.3.15) 

These equations show that \zi{t), — >is a decaying state and ,in fact all states 
oi 4>^\1-L are decaying states, while \zi{t),+ > is a growing state, and all 
states oi (j)^\Ti. are growing states. 
It can be proved that [15]: 

K\zi - >= 1^1,+ > 
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K\zi,^ >= \zi,- > (5.3.16) 

which is a natural fact, since growing states must be transformed in decaying 
states by the time-inversion operator and vice versa. Therefore we have: 

K-.r-^ <t>i 

K : ^ (j)"! (5.3.17) 

The following limits are valid (cf. eq. (5.3.12") : 

lim l^i, - >= lim \zi, + >= |1 > (5.3.18) 

A— >oo A— >oo 

Therefore — > and > can be considered as version of the unstable 
state |1 > in spaces 0^ and 0^. In fact, the difference between these vectors 
and |1 > is a 0(A), since when A = the interaction disappears. Let us 
remember that the survival probability of state |l(i) > was: 

P{t) = I < l|l(t) > |2 =< l|l(t) >< l(t)|l > ....(5.3.19) 

P{t) shows the initial Zeno effect behavior, then an exponential behavior and 
finally the oscillatory Khalfin effect behavior. If we make the substitution 
|l(i) >^ l^i(i), — > we obtain: 

P{t) P'{t) =< l\zr{t), - >< zi{t), -|1 >= e-^i* (5.3.20) 

and only the exponential behavior remains. Thus the physical nature of the 
state 1^1, — > would be the one of a decaying unstable ideal state, where 
we have ehminated the Zeno and Khalfin effects, because these effects are 
contained in the last term of the r.h.s. of eq. (5.3.12") (also called "the 
background"). Namely, the three effects are mixed if we use the time evolu- 
tion based in eq. (5.3.6), but Zeno and Khalfin effects can be separate, from 
the exponential behavior, if we use the evolution based in eq. (5.3.12). Eq. 
(5.3.20) show also that 7jf ^ is the mean fife time of the unstable states. Fine 
graining can be thought as an approximation of real states that eliminates the 
unimportant Zeno and Khalfin effects. Zeno effects is unimportant because 
it takes place at t = while, we are generally interested in the phenomena at 
t — > oo. Khalfin effect it is uninteresting because essentially it is an oscilla- 
tory effect, around the exponential behavior, while we are interested in mean 
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values only. But as we have said in the introduction it is not completely clear 
if these ideal exponential states are just, mathematical useful,effective state 
or real physical states. 

Furthermore, fine graining school need to work with rigged Hilbert spaces 
0^ and 0^ to solve the problem of the arrow of time as we shall see in the 
next section. 

Friedrichs model is just an example, but its rigged Hilbert space structure 
can be found in every scattering process [38] . Therefore even if we will base 
our reasoning in this model, what we will explain below is rather general. 

5.3.4 Mixed states. 

Let us now introduce the arguments of the next subsection written the evo- 
lution equations of mixed states in our model. A mixed arbitrary state at 
time t = can be expanded in basis {|a;, ret >} as: 



We can as well use the advanced basis, but this is all what we can say in space 
H. But in space 0^ we can use the basis >,\u!,— >} (introduced in 

eq. (5.3.12')) and expand p as: 




(5.3.21) 



and its time evolution reads. 




(5.3.22) 



P = Pll\zi, - >< Zi, -| + J^{piu;\zi, 

+Pcvi\(^, - >< zi, -\)du;+ 



Xuj - 



...(5.3.23) 



and its time evolution reads: 



P{t) 



p*(i) + e-5^iVi(i) + e-^^V2(i) 



(5.3.24) 



where: 




(5.3.25) 
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and: 




Now since 71 > 0, pi(t) oscillates, and, P2{t) is time invariant we have: 



which seams very close to the strong limit we are looking for. The only prob- 
lem is that it is not a dynamical equilibrium state, since it oscillates. 
Nevertheless from the thermodynamical point of view is a thermody- 
namical equilibrium state, since it has a constant (and maximum) Gibbs 
entropy. In fact , it is evident (from the quantum version of theorem (3.6.2)) 
that Gibbs entropy is constant for the time evolution (5.3.22), therefore it is 
constant for the time evolution (5.3.25) which it is similar. It is clear that 
is all what we can ask to the model, since the field cannot go to dynamical 
equilibrium because the modes of the field are uncoupled. 

Therefore, from the thermodynamical point of view eq. (5.3,28) reads. 



and it is the strong limit we are looking for. As in the case of eq. (4.5.1.11) 
this limit belongs to the corresponding rigged Hilbert space. 

Where is the miracle that allows to pas from the oscillatory evolution(5.3.22), 
with no limit to the partially dumped evolution (5.3.24) with a thermody- 
namical limit? The miracle is that eq. (5.3.22) is valid in space C = Ti xH. 
while eq. (5.3.24) is valid in space $^ = q!)^ x 0^ so really eq. (5.3.28) is a 
functional equation that can be interpreted as: 



where 0_ is an operator of the test operator space $_ = 0_ x the space of 
measurement operators we have chosen. Therefore the miracle happens just 
because we have chosen a convenient test space for our physical measurement 
apparatuses. Eq. (5.3.30) is a weak limit that is similar to the weak limit of 
mixing classical states and a consequence of theorem 5.2.1 since our model 



lim {p{t)-p,{t))^Q 



(5.3.28) 



hm pit) — p* 

t — ^00 



(5.3.29) 



hm (p(i)|0_) = (p,|0_) 



(5.3.30) 
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has a continuous spectrum. We will continue this line of reasonings once the 
more complete example of the next subsection would be introduced. 

In the practical case we will study in section 8 the field \cu > will be the 
thermic radiation field within the universe that we can considered as thermal- 
izcd., from its beginning, by other interactions than those of eq. (5.3,4), then 
it can be classically chosen as a Boltzmann thermic distribution function: 

p^^ZT'h'^ (5.3.31) 

where T is the temperature and Z a normalization function, and the dumping 
terms are produced by nuclear reaction phenomena within the stars and 7f ^ 
is the characteristic time of these nuclear reactions. We will use this model 
in section 8. 

5.4 Friedrichs model for many oscillators. 

We will now introduce a not very realistic physical model, that nevertheless, 
is the simplest one for our purpose. Let us consider a set an infinite (or a great 
number) of uncoupled harmonic oscillator, labelled by uj , with hamiltonian: 

H^ = Lj{ala^ + ^) (5.4.1) 

where aj^ and a^^ are the creation and annihilation operators of the harmonic 
oscillator. The total hamiltonian reads: 

POO 

H= / H^duj (5.4.2) 

Jo 

H,^ can also be written: 

^^c-E^ir^ (5-4.3) 

n 

where: 

H^"-^ = u{n\n,u} >< n,u}\ + -) (5.4.4) 

|n, cj > is the a;-oscillator in the n exited state (n = corresponds to the 
ground state) 
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Let us suppose that each of these states is coupled with a field represented 
by a set of infinite states \n,cu,w > in such a way that now the coupled Hj^J^'^ 
reads: 

ifj^") = Lij{n\n,uj >< n,u\ + |) + dw\n,uj,'w >< n,uj,'w\ + 
A /q°° dwgn,u){w){\n, cu >< n,u!,w\ + \n, cu, w >< n,u!\) (5.4.5) 

where A is a coupling constant and gn,uj{w) a coupling function which neces- 
sarily has the property gQ^uj{w) — 0, since the ground state of each oscillator is 
stable and therefore it is not coupled with the corresponding field that would 
produce its instability. Therefore we have constructed a model which can be 
considered as a infinite repetition of the Friedrichs model of last subsection. 
In this non-realistic model the instability of all the states, with the exception 
of the ground states, is obtained by couphng a field to each oscillation mode. 
It is, in fact, a non-economical way, but it serve to our purpose, which it is 
now only to find the laws of unstable evolutions. 

Now, using the procedure of the previous subsection in each Friedrichs 
model of each mode we can diagonahze each operator Hj^^ and we will obtain: 

i/W ^ 4'^)n|n,a;,+ >< - | + ^o; (5.4.6) 

where, for simplicity, we have omitted the field term, where z';^^ — uj, since the 
ground states are not perturbed, and, Imz^^ < forn 7^ 0. If we renormalize 
and eliminate the |-terms we obtain: 

/•oo 

/ y^z^''\\n,uj,+ ><n,uj,-\duj (5.4.7) 

where we have put the factor n inside z^^ . 

Let us now consider a density matrix p — p(0) e C — Ti. x H, that can 
be expanded in basis {\n,u!,— >} as: 

POO 

P= y2p"',n';io\n,u;,- ><n,u;,-\du; (5.4.8) 

We will always work with these density matrices below. And these is the 
essential fact. Since \n,uj,— >G 0^ and, therefore p G x 0^, what we 
have done, in choosing the expansion (5.4.8), is to assume that our operator 
space is = 0_ x 0_ in such a way that now we have the Gel'fand triplet: 

C £ C (5.4.9) 
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and, therefore, p G $^ = 4>^ x 0^. On physical ground what we are doing is 
to postulate that our measurement apparatuses correspond to operators in 
We will discuss this postulate below, but we can see immediately that 
this is the price to pay to get the strong limit we are looking for and the 
corresponding unstable time evolution. 

In fact, from eq. (5.3.12) we can obtain the time evolution: 

p{t) = e-^^V(O) = 
CEnPn,n';.,e-^(^^"'-^^"''*)*|n,a;,-><n',a;,-|da; (5.4.10) 

but since 

-z(4") - z^f>) = - i3(f^) - '-{^1^^ - 7(7')) (5.4.11) 

and 7i"\7<lr'^ > 0, (only = 0), thus: 

/•oo 

lim p(t) = / ^po,o;c^|0,c;,- >< 0,t,;,-| = (5.4.12) 

and we have obtained our strong limit, equivalent to (5.2.9). Furthermore, 
now we have the time evolution to this limit, eq. (5.4.10). To obtain this 

result we have used an infinite set of continuous field that we neglect in all the 
formulas above. Some how we have "traced away" these fields. But the result 
will not change, from the physical point or view, if we write all these fields. 
The result we have obtaining regarding the states of the harmonic oscillator 
will be the same, these oscillators reach to the equilibrium showed in the last 
equation, but the fields will continue to oscillate and they will be always far 
from the equilibrium (like in the last part of the last subsection). This is not 
surprising, since these fields have no self-interaction or mutual interaction, 
therefore they cannot reach to equilibrium. Then, to neglect these field was 
only a useful shorthand with no physical consequences- (provided we take into 
account all the warnings we made in the last part of the last subsection). We 
must also remember that the quantities of eq. (5.4.12) are just functionals 
over the space $_ thus if we contract this equation with any vector of this 
space we will find the weak version of limit (5.4.12) showing that in this 
example theorem 5.2.1 is fulfilled. 

If we collectively call 27 to all the 7's or 27 is the inverse of the char- 
acteristic life time of the system or if we call 27 to the smaller of them to 
maintain the leading term only, eq. (5.4.11) reads: 
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= + pie-^* (5.4.13) 

as usual we have: 

trp = trp^ = 1 (5.4.14) 

.since matrix p is the usual one and matrix is an expansion of stable states 
(5.4.12) and norm must be conserved (cf. Appendix 5.B). But: 

trpi = (5.4.15) 

as a consequence of eq. (5.3.14'), showing that this matrix is something like 
a fluctuation around the equilibrium state. 

Let us finally observe that in this model we can not pretend that 
would be the equilibrium state of the canonical ensemble. To obtain that 
result obviously we must coupled the oscillators-among themselves and the 
model will be much more complicated. To mimic a canonical ensemble at 
temperature T in this model the best we can do is to make the following 
choice: 

Po,o;u> = (5.4.16) 

with this choice we have the correct equilibrium distribution and the field 
produce the irreversible evolution toward this equilibrium. The evolution of 
p{t) then reads: 



PC 

Pit) = / 
Jo 



^e-^^p^f + e-2-V(^)pS'^^ 

1 2 



du (5.4.17) 



where /(a;) -is an arbitrary function, so the conclusions are essentially the 
same as in the last subsection. 



5.5 Appendix 5.A.-Hardy class functions [29]. 

A complex function /(cu) on it!+ is a Hardy class function from above (below) 
if: 

1.-/(0;) is the boundary value of a function f{z) of the complex variable 
z = X + iy that is analytic in the half plane y > (y < 0). 

ii.- j'^^ \ f{x -\- iy\^dx < k < 00, for all y with < 7/ < 00 (—00 < y < 0). 
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5.6 Appendix 5.B. Computation and conservation of 
norm, trace, and energy. 

The trace of a density matrix p is: 



and it is invariant under changes of basis. Making the same procedure that 
we used to go from eq (5.3.6) to eq. (5.3.12) we can obtain: 



As the fine graining theory deals with states that vanishes when t — > oo, we 
can be worried guessing if the norms, the traces, or the energy are conserved 
in this theory. There is no problem since we can state the following results. 

i. -Prom eqs. (5.3.14,14') we see that unstable density matrices hke \zi, — >< 
— I or 1^1, + >< 2^1, +1 have null trace this fact is possible since we are not 

working in Hilbert space.. 

ii. -Using eqs. (5.3.15) and (5.3.14') we can see that the trace of eq. (5.B.2) 
is conserved as the usual trace of eq. (5.B.1) 

iii. -If the trace is conserved also the norm of pure states is conserved. 

iv. -The mean value of the energy in one of these unstable states hke 
\zi, — >< Zi, — I reads: 

< H >=< zi, -\H\zi, - >= zi < zi, -\zi, - >= (5.5.3) 

and therefore the energy of the states that vanish when i ^ oo it is zero, 
creating no problems with energy conservation.. 

6 Coarse- Graining and Trace. Time Asym- 
metry. 

6.1 Coarse graining. 

Let us go back, for a while, to the classical regime. Usually coarse-graining 
is based in the fact that the dynamical variables cannot be measured with 
infinite precision, i.e., there is always an error and also we cannot compute 




(5.R1) 




(5.S.2) 
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with an infinite number of digits. Perliaps tliere is a fundamental graininess in 
nature but tliis graininess it is not yet neitlier tlieoretically nor experimentally 
found. 

Coarse graining can be introduced by partitioning the space X into finite 
(or discrete) number of cells Ai that satisfy: 

[jA^x, A n = 0, .if..i ^ j (6.1.1) 



This partition is arbitrary but it must be non-trivial with respect to some 
measure /i namely: 

< pL{Ai) < pl{X) (6.1.2) 

for all values of i. For every density p within each cell Ai of the partition, we 
can compute the average of p as: 



< P >i 



1 



p{x)ii{dx) (6.1.3) 



p{Ai) J A, ' 

and the coarse grained density respect to the partition is given by: 

P(x)=J:<p>, U(x) = {Y:-Jj-\1a,){U,\}\p) = np(x) (6.1.4) 

where 1^. is the characteristic function of the cell A^ and 11 is the projector 
defined by the partition. 11 is a projector since: 



E 



1 



E 



A*(A) 



|1aJ(i. 



E 



|lA,)(U, 



^ lU)(Uj//(A)<^i,- = Y.^^\^a.){Ia,\ = n 



From the reasonings of the introduction or from theorem 4.3.1 we can deduce, 
in the case of finite number partition (for discrete number see [7] ) the 

Theorem 6. 1.1. -If Pt is a p*-mixing Markov operator with a unique 
stationary density and {A} is a non-trivial partition of the phase space 
X, then: 



lim iPtpY =limP* (6.1.5) 

t— +0O t— >oo 
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for all initial densities. 

Thus we have obtained our coarse graining strong limit. Now we can 
consider the transformation: 



from eq. (6.1.5) we see that this transformation has a strong limit (i.e. in 

the norm) and therefore, according to theorem 4.4.1 it is exact. 
Then using theorem 4.4.2-we can say about entropy that: 
Theorem 6. 1.2. -If Pj is a reversible mixing Markov operator with a 

unique stationary density and {Aj,} is a non trivial partition of the phase 

space X, then: . 



for all initial densities p. 

But we must realize that the way the conditional entropy converge to zero 
depends on the way in which the coarse graining is carried out [6]. It can be 
proved that the rate of convergency of entropy to equilibrium becomes slower 
as the measurement technics improve and the coarse graining becomes finer 
(!!!). Such phenomena have not been observed. Thus is most unlikely that 
trivial coarse graining would play a role in determining thermodynamical 
behavior, if a natural graininess it is not found. Candidates for this natural 
and universal graininess would be: 

i. - The graininess produced by operators n, 11 which are introduced using 
fine graining methods [15], [40]. But, really, this is only the coarse-graining 
old version of the fine-graining method. 

ii. - The graininess introduced by the universe event horizon [41]. 

iii. -The graininess introduced by Plank's quantities, e.g. it seems that it 
is absolutely impossible to measure length smaller than Plank's length. 

But the physic related with these graininess is still under research. 

More general projectors than those defined in eq. (6.1.4) can be used, as 
we have seen in the introduction, since any projector will do the job done 
in eqs. (1.2.1.2) and (1.2.1.3). A theory that uses one of these generalized 
projectors will be called, by extension, also a "coarse-graining" theory. 

Coarse-graining can be used also in the quantum case. Then 11 is a 
projector over the quantum Liouville space C = TC <^TC. Using the quantum 
theorem (5.2.1) we can obtain the same results as in the classical case, where 
we use the theorem (4.3.1). 




(6.1.6) 



Jim Hcmpit)r\p*) = o 

t— »oo 



(6.1.7) 
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6.2 Time asymmetry in coarse graining theories. 

All the limit of the previous section were computed as t +00 but it is clear 
that all these limits are also valid when t —00. Therefore in coarse graining 
theories there is equilibrium both in the far past and the far future. This 
fact can be easily verifier with the baker transformation, where we have a 
set of infinite lines, horizontal for the far future and vertical for the far past, 
that will be considered as a uniform equilibrium distribution function for 
any coarse graining partition. It is also evident that if the initial distribution 
function has an adequate symmetry, e.g. the one of the characteristic function 
of a square domain, the evolutions, toward the past and toward the future will 
be strictly symmetric, but this will not be the case if the initial distribution 
function is not symmetric (other calculations about the baker transformation 
behavior can be found in ref. [6]) 

Going now to the quantum case we see a quite similar phenomena, even 
with no coarse graining. Let us consider the state |1 > of Priedrichs model 
(which can be considered as a symmetric initial condition, as the character- 
istic function of an square domain, in the case of the baker transformation). 
The behavior of the survival probability P(t) , as shown in fig.2, is completely 
symmetric with respect to t = 0. Thus classically if we use coarse graining 
technics, or quantum mechanically if we use only states of the Hilbert space 
we will find that past is only conventionally different than future. What is, 
then, the way to distinguish past from future? It is the method that we 
explained in the introduction: Take the time t = 0. Consider the set of 
evolutions of the system for t > (for all possible initial condition) and let 
as call it H- in the quantum case (or £_ in the classical case) . It is identical 
to the set of evolutions for t < (for all possible- initial conditions) , that we 
shall call H+(or £+). The existence T or K, the mathematical transforma- 
tion, that relates the future evolutions with the past evolutions shows that 
these sets of evolutions are identical. In fact. 

Theorem 6.2.1. For every evolution p{t) G C_ t {t > 0) there is a 
time symmetric evolution p{—t) e {t < 0) if the evolution equation are 
reversible. 

Proof: From the definition of a reversible evolution (2.3.21) for every 
p{t) G jC- there is a physical evolution p{—t) G defined by 

p{-t) ^lCp{t),....t > 0, □ 
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Also we have: 

T -.C^^ C+, T : £+ ^ £_ 

These two sets of evolutions £_ £+(or 7i-H+) are the two mathematical 
structures that we introduced in subsection 1.1. As they are identical (cf. 
eq. (6.2.2)it is irrelevant to chose one or the other. So let us chose one of these 
structures to build our theory, let us say and forget the other. Now we can 
say that the theory begins at t = and goes toward the future for t > (or 
toward the past since the choice of one word or the other as just conventional 
as the choice between £_ and £+). It is quite evident that this theory 
developed in the lapse < t < oo will fulfill all our requirements, provided 
we forget all about the lapse — oo < t < 0. These are the characteristics of the 
resulting theory if we use coarse graining and usual Hilbert space quantum 
states. Even if successful in many respect it produces, should we say, a 
certain uneasiness. 

6.3 Traces. 

Let us now consider the classical case only. 

Let X and Y be to topological Hausdorff phase spaces, ip : Y ^ X a given 
continuous function on X , and St : Y ^ Y a dynamical system operating 
in phase space Y. A function h : R ^ X is a trace of the dynamical system 
if there is a point y in space Y such that h(t) = Lp{St{y)), for all times t (this 
meaning of the word trace must not be confuse with the one we use when we 
speak about the trace of a matrix) . 

It can be proved that every continuous function in a space X is the trace 
of a single dynamical system operating in a phase space Y, therefore we have 
the quite surprising 

Theorem 6.3.1.-[6] Let the phase space X be an arbitrary but topolog- 
ical Hausdorff space. Then there is a second phase space Y also topological 
and Hausdorff, a dynamical system St operating in Y and, a continuous func- 
tion ip : Y ^ X such that every continuous function h : R ^ X is the trace 
of St- (A topological space is Hausdorff (or separable) if any two distinct 
points possess disjoint neighborhoods) 

That is, for every h there is a point y in phase space Y such that h{t) — 
(p{St{y)), for all times t. 

Let us now consider the trajectories of a dynamical system: If we have a 
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dynamical system St operating in a phase space Y , then only three possible 
types of trajectories may be observed: 

i. - The trajectory is a fixed point such that SfX^ — x^, for all t. 

ii. -The trajectory is a non intersecting curve, with the property St{x) ^ 
Stix) ift^t'. 

iii. -A periodic trajectory such that St{x) = Sf+rix), for all times t, being 
T the period. 

But nothing prevents the existence or non periodic intersecting trajecto- 
ries h(t), in space X if ip : Y ^ X. Thus we can demonstrate the following 

Theorem 6. 3. 2. -Let the phase spaces X and Y be topological Hausdorff 
spaces and h : R X an intersecting and non periodic trace of a dynamical 
system St : Y ^ Y . Then the entropy of densities evolving under the action 
of h is either constant or increasing. 

Proof. The proof is based on the trivial observation that if h is intersecting 
and non periodic, then at every intersection point x on the trajectory h 
the inverse h~^{x) is not unique. Therefore the trace h is the trajectory of 
a semidynamical system, and since semidynamical systems are irreversible, 
from theorems 3.6.1 and 3.6.2 the entropy is either constant or increasing. □ 

Thus the simple act of taking a trace of a dynamical system (with time- 
constant entropy) may be sufficient to generate a system in which the entropy 
is increasing. But for certain class of traces much more can be said. 

Let X and Y be two different phase spaces with normalizable measures 
//* and i/* and associated densities p* and cr*, respectively, and Tt : X ^ X 
and Sf : Y ^ Y he two measure preserving transformations. If there is a 
transformation (f :Y ^ X that is also measuring preserving, i. e. if 

iy,{if-\A)) ^ li^A) (6.3.1) 

for all subsets A of the phase space Y, and such that o 93 = (po St, then Tj 
is called a factor of St- Prom this definition the trajectory of the factor Tt is 
a trace of the system St- Then we have the following 

Theorem 6. 3. 3. -[42]. Every p*-exact transformation is the factor of a 
Kolmogorov automorphism. 

This theorem precise the things we must do if we want to find an exact 
transformation with all its nice properties: 

i.-We must show that the system we are working with is a Kolmogorov 
system. This can be difficult from the mathematical point of view, but as 
chaos is very frequent in nature it is not a very restricting physical condition. 
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ii.-Then, according to theorem 6.3.3., every measure preserving factor 
will produce an exact transformation. The problem is just to find the most 
convenient one. 

As an example let us consider again the baker transformation. It can 

be prove that this transformation is a Kolmogorov automorphism, endowed 
with a constant entropy. However the system corresponding to coordinate x 
is a factor of baker transformation. Also it is identical to the dyadic Renyi 
transformation: 

T{x) = 2x..{modl) (6.3.2) 

which is uniformly exact and whose entropy smoothly increases to zero by 
theorem 4.4.2 

We have show that coarse graining produces no substantial difference 
between past and future. This is not the case with traces, as we can see from 
baker transformation where that x-side of a parallelogram will always increase 
toward the future and decrease toward the past, Thus coarse-graining do not 
produces time asymmetry while traces do produce this phenomenon. 

Let us now give all the whole panorama : 

-We have projectors 11 like those introduced in section 6.1, namely such 
that: 

n : £y ^ Cx, = n, (6.3.3) 

where Cy — JC is the state space and Cx is the space of relevant states. 11 
has not inverse H"^ since n^n~^ = nil""^ yields H = 1. 
-We have traces: 

ip-.Y ^X, (6.3.4) 

namely mapping between phase spaces, (p can have an inverse, and in this 
case (p{Y) is dense in X, or it do not have an inverse when the dimension of 
X is smaller than the dimension of Y, like in the case of eq. (6.3.2) 

Finally let us remark that traces define a mapping in the corresponding 
Liouvillc spaces. Let Cx and Cy be the corresponding Liouville spaces to 
the phase spaces X and Y. Then to the mapping ip : Y ^ X corresponds 
the mapping: 

A-^ : Cy l^x (6.3.5) 

(the -1 is just a matter of convention) defined by: 

A"V(|/)=P(^"'(^)) (6.3.6) 
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In the next subsection we will study even more general mappings. 
6.4 Generalized traces. 

In section 6.3 we are forced to work in the classical case only, since we have 
used phase space. Now we would like to generalize the trace notion in order 
to work also in the quantum case. 

A generahzed trace is given by eq. (6.3.3) if eq. (6.3.4) is not fulfilled, 
i. c. it is a mapping between Liouville spaces not originated by a mapping 
between phase spaces. Being a mapping like (6.3.3) it is something like a 
"projector with an inverse". But now spaces Ly and Lx can be classical or 
quantum Liouville spaces. These generahzed traces are typical of the fine- 
graining formalism, and try to show it as a kind of generalization of the 
coarse-graining one. 

Let us consider the particular case Ly=^^^ Lx = L. Let us also consider 
the basis {|1 >, jo; >} of eq. (5.3.1) that we shall call {|i >} and such that 
H\i >= Zi\i >,Zi e C. let us define the basis {|ii)}, = \i >< i|-Let us 
also consider the basis {l^i, — >, \u;.— >} of eq. (5.3.23), that we shall call 
>} and in the same fashion let us define the basis {\ij, —) — 

>< Using the basis >, + >} we can, as well define a 

basis {\ij, +)}. Then we can define a generalized trace as: 

A-^ : ^ jC (6.4.1) 

A = E -'A-' = E +1 (6-4.2) 

ij ij 

Namely A is the transformation that make correspond each state p of 
space £ to a functional in space A looks like just a "change of basis". 
But really A is much more than a change of coordinates since it takes vectors 
of one space to vectors in another space. Therefore to weak limits in C 
corresponds strong limits in $^ and the generalized trace A embodied the 
solution of our problem: to go from weak limits to strong limits and can be 
considered as the symbol that synthesize the fine graining technic. 

Some observations are in order: 

i.- Since A is a generalized trace, therefore as a trace it contains time 
asymmetry. In fact A defined in eq. (6.4.1) is related with dumping phenom- 
ena that produces equilibrium toward the future and should be called A_. 
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We well, define a A+ related with creation phenomena that implies 

equilibrium in the far past, namely: 

A+ : ^ C (6.4.3) 

A+ = ^ \tj, +){tj\, A;^ = ^ \tj){tj, +1, (6.4.4) 

ij ij 

if we chose A_rather than A+ we are creating a time-asymmetry. 

In order to see the relation of the two A's let us introduce the star- 
conjugation: 

A* = 1CA^IC^ (6.4.5) 

then it is easy to see that: 

A+ = A*,....A_ = A* (6.4.6) 

a_a:^ = a+a; = i, a:^ = a* = a+, a;^ = a; = a_,... (6.4.7) 

for these last equation we can say that the A's are star-unitary. 

ii. - Using generalized trace A we do not lose any information, (in the case 
of usual traces we lose information if dimX < dimY, as in the example of 
the baker transformation before eq (6.3.2) but in the case of the A-trace the 
dimension of the two spaces is the same and C is dense in $^). But it can 
be demonstrated that this generalized trace A, some how, renormalize the 
infinite amount of information contained in £ [43] 

iii. -A generalized trace is not a trace, so there is not mapping in the 
corresponding phase spaces that make trajectories correspond to trajectories. 
In this sense using fine-graining technics trajectories loose al their importance 
and even have no meaning. [44], [45]. 

iv. A-trace allows to define a Hilbert space where the time-evolution are 
irreversible.. 

In fact: using the bases we have introduce we can deduce that: 
l = El^i' -)(U,+I (6.4.8) 

ij 

L^J2i^,-z;M-){ij,+\ (6.4.9) 

ij 



77 



where as ImZi = — y < 0, it is Im{zi — 2;*) < 0, The time evolution operator 
is U{t) = e~*^*and UW = 1, i. e. U is unitary. Let us now define a modify 
huoviUian: 

G = A*LA = J2i^, - z*)\ij){ij\....{6A.10) 

which induce a evolution W{t) = e"*"-^* such that WW"^ 7^ 1, and, therefore 
it is not unitary but star- unitary WW* = 1. the two evolution are related 
by: 

W{t) = A*U{t)A (6.4.11) 

We can also define A-density matrices as, related by the A-trace (6.4.1), as: 

PAit) = A>(t), ...p{t) = ApA(t).. (6.4.12) 

where p{t) e $^,pA(i) e C that evolve as: 

p{t) = U{t)p{0), ..pA{t) = W{t)pA{0) (6.4.13) 

Eq. (5.3.24), translated to the pa language reads: 

PA{t) = PA. + e-^^^VAi(t) + e-^''pA2{t), (6.4.13') 

that can be also obtain from eq. (6.4.13), since the p\{t) evolve under the 
action of the operator e~^^^ and G has complex eigenvalues (cf. (6.4.11). 

Then the space of the p's can be considered as an ideal reversible world 
of reversible equations, namely the ideal world of Newton, endowed with 
unitary evolutions, (David Bohm would say that this is the space of implicate 
order [43]), while the space C of the p^ is the real, physical, irreversible world 
of Boltzmann, endowed with non-unitary evolution (just star-unitary, David 
Bohm would say that this is the space of explicate order [43]). Between 
these two worlds A establishes a canonical mapping (David Bohm would say 
a "metamorphosis" [43]). Even if the p^s live in the ordinary Liouville space 
they evolve with a non-unitary law (cf. (6.4.13)), so A— trace achieve the 
dream of physicist: it creates an ordinary Hilbert space where the evolutions 
are non-unitary and irreversible. Precisely, the essence of the fine-graining 
formalism was to maintain the time symmetric primitive equations (with 
operator evolution U (t)) and to obtain time-asymmetry by choosing a typical 
time-asymmetric space Cy = A-trace change these roles We get a time- 
asymmetric equation (with evolution W{t)) in a time-symmetric space C as 
in the coarse-graining case. But, of course, the physics remain the same. 
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v.- Using the A Lyapunov variables can be find very easy since: 

{p{tMt)) = (p(0)|C/tC/|p(0)) = (p(O)Ip(O)) = const, (6.4.14) 

therefore it is not a Lyapunov variable, but: 

Y = {pA{t)\pA{t)) = {pj,{t)\WW\pj,{0)) = var., (6.4.15) 

or 

Y = (p(t)|(A*)tA*|p(t)) = {p{t)\M\p{t)) = var., ...M = (A*)tA*, ...(6.4.15') 
Precisely: if 

Ipa(O)) = ^PiiKj) 
the corresponding time evolution is: 

ij 

but Im{zi — Zj) — —Tij < 0, so: 

{pA{t)\pA{t)) = Y.\pij\'^~''''' (6-4.18) 

ij 

is always decreasing and it is, therefore, a Lyapunov variable. 

In more general cases, than the one of cq. (6.4.1)-(6.4.2), it can be proved 
that every rigging corresponds to a A-trace and vice versa [43] 

6.5 Time asymmetry in fine graining theories. 

Let us begin computing the conditional entropy Hc{p\p*) (cf. eq. (3.6.1,2)) 
in the case of the classical evolution (4.5.1.12) (in the quantum case we 
have the time evolution (5.4.13)). If we want to use the classical equation 
for He and we have a quantum density matrix we must first transform this 
quantum matrix to the corresponding classical distribution function, using 
eq. (6.A.1). But we can use directly the definition of He if we define the 
logarithm of a quantum density matrix as the operator whose eigenvalues 
are equal to the logarithms of the eigenvalues of the primitive operator [3]. 



(6.4.16) 
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As Wigner integral is linear the quantum analog of eq. (5.4.13) is the same 
equation). Then: 



Hc{p\p*) = - / (p* + e-^Vi) log(l + e-^'^)dx (6.5.1) 




considering that \pi\ << or t » 7 ^ we can expand the logarithm and 
since trpi = 0, taking into account eq. (6. A. 3) we have: 



which is negative, growing, and, with a vanishing limit when i — > 00, so it 
accomplish all the properties to formulate the second law of thermodynamics 
in its third order form. But we have obtained this satisfactory conclusion 
because we have work with the operator test function space $_ and the 
quantum states belong to space $^ (albeit some mathematical problems, 
since we are computing the log of a vector of a rigged Hilbert space; these 
problems can be solved, in principle, if we use the generalized trace A, of the 
previous section, and if we substitute the p's by p^s, since these last density 
matrices belong to C ,so they can be used with no problem, but they keep the 
evolution properties of the p's, namely the dumping factors of eqs. (6.5.1) 
and (6.5.2) as in eq. (6.4.13'); namely we define 



. So we can go now to the central problem of the origin of time asymmetry 
in fine graining theories. 

Let us now consider an isolated system which is all our universe, there 
is nothing we can know about the exterior of the system and the system 
cannot interact with something out side the system. If the time evolution 
equations of a theory are time symmetric it is quite impossible to brake this 
symmetry by rigorous mathematical manipulations, symmetry will always 
appear, some way or another. Nevertheless, the examples we gave show that 
normally in these theories we can find two extensions of Liouville space C 
They are the rigged Hilbert spaces and $^ which are defining using the 
test spaces $_ and (usually these test spaces are nuclear spaces that can 
be consider the spaces corresponding to the operators of the measurement 




(6.5.2) 





Pa log — dx) 

Pa* 
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devices, as explained in appendix 4. A). Time symmetry makes that these 
spaces-are related by: 

/C ^ /C (6.5.3) 

and, therefore, they are identical. To chose one of the other is an irrelevant 
choice. As irrelevant as to trough a dice with the same number in all its 
faces. So if we chose one space or the other, physics do not change. Both 
spaces are only conventionally different. Any possible difference could came 
only for the exterior of the system and there is nothing there to interact 
with. Nevertheless in each space or future is substantially different 
than past, since there is equilibrium toward only one of these directions and 
we can call this direction the future. So let us chose one of the spaces, we 
then establish a time asymmetry and we can formulate the second law of 
thermodynamics, as we have done, and our problem is solved (compare the 
solution with the coarse graining case, it is not so different) 

We can say the same talking about the choice of the generahzed traces 
A_ or A+ and work within the Liouville space C, the space of physical states. 

The same trick can be done in various different ways e. g.: 

i. -In papers [26], [27], [28], and [32] two semigroups are defined, each related 
to a rigged Hilbert space, and one of these semigroups is arbitrarily chosen. 
One semigroup is obtained expanding the solution of the evolution equation 
in a basis of and the evolution turns out to be well defined for t e 
(— oo, -|-cxo], namely it is not well defined for i ^ oo. The other semigroup 
has the inverse properties. 

ii. -In the book [40] a projector n_ = J2i l^h +1 is defined and con- 
sider as the projector on the really relevant space. But n+ —KH is identical 
to n_ , so we must chose one or the other as in the previous cases. 

So in all these cases we must do a conventional choice to find a mathe- 
matical structure: a space, a semigroup, a projector,... such that using this 
structure the future exhibit substantially different properties than the past. 

Someone may say that we have not explain time asymmetry, since we have 
just introduce it by an arbitrary choice. To answer this criticism we must 
remember that physics really never explains. It merely find the mathematical 
structure more adequate to foresee the physical phenomena: c. g.: the more 
adequate mathematical space, the more adequate mathematical equations, 
etc. The curvature of space-time do not explains gravity, it happens that a 
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riemannian manifold is the best mathematical structure to deal with gravity. 
Analogously, it simply turns out that the most adequate mathematical space 
to explain time asymmetry and the second law of thermodynamics is a rigged 
Hilbcr space not the usual Hilbcrt space: 5*0 the relevant important choice 
is between the mathematical structures C or $^ (or which is the same think 
$!^) The choice between these last two rigged spaces, or is on the 
contrary, irrelevant and physically unimportant 

6.6 Comparison between Fine-Graining and Coarse- 
Graining. 

As we can see coarse-graining and fine-graining are very similar. 

Both are based in theorem 4.3.1 about the weak limit of mixing evolutions. 
Coarse-graining obtains a strong limit via a projection, fine-graining obtains 
the "strong limit" using functionals. 

Both obtains their arrow of time defining a pair of time. symmetric struc- 
tures. The pair of t > 0, and t < evolutions in the case of coarse- 
graining. The pair of rigged Hilbert spaces in the case of fine-graining. 
In both methods one of these structures is conventionally chosen. 

The main weakness of coarse-graining is that the projector is not defined 
in a canonical way. 

The main weakness of fine-graining is that we are force to enlarge the 
space and we do not know the exact nature of the objects we must add. Are 
these ideal unstable states just mathematical usefuU tools (like Fadeev-Popov 
ghost) or real physical objects? The answer to this question depend on the 
point of view that we would take studying the problem. In fact: 

i. Any decaying state was always created by a creation process. The 
quantum state that corresponds to the creation process followed by the de- 
caying belongs to H (like the vector |1 > of the Priedrichs model with the 
survival probability (5.3.19), the one of fig. 2). Nevertheless, if the hfetime 
of the decaying state is very large, we use to neglect the creation process and 
to consider the state just like a decaying state with exponential decaying sur- 
vival probability (as in eq. (5.3.20)). This is the state \zi, — > that belongs 
to 0^. So the quantum theory, that uses these states, could be considered 
as an effective theory, where creation process are neglected. We can say the 
same for classical theories. In Bakers transformation a regular density, with 
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a regular support, will have a creation process and a symmetric decaying 
process, towards equilibrium, much in the same way as state |1 > • But if 
we study the time evolution of a "horizontal Dirac comb" states we will find 
that these ideal states have no creation process, as the state \zi,— > . 

ii. Nevertheless, what could be just an usefull simplification, when using 
states with large lifetime, can be a rigorous fact in the case of the universe, 
where, as we do not know its creation process, this process must be necessarily 
neglected. 

iii. So the new unstable states added to physical space are similar to 
planes waves, they are eternal objects with no creation process and, in fact, 
if we define plane waves in a rigorous way we need a rigged Hilbert space 
to do it. From this point of view, coarse-graining physicist would be like 
stubborn persons that only work with waves packet and refuse to use plane 
waves because they "are not physical objects". 

iv. If we allows time go to infinity and we would like to consider the 
rigorous equilibrium state at infinite time, this state belongs to $^ as in the 
case of Baker's transformation, so we are forced to work with a fine-graining 
theory. But if we content ourselves with approximate equilibrium states at 
finite time, arguing that i — > oo is physically impossible, we do not need 
these states. 

So the real nature of the new states is open to discussion. However a 
fine-graining physicist can take a conservative attitude and consider the new 
states just as ideal states, namely just as useful mathematical devices as 
plane waves are. 

Is it the choice of the fine-graining or the coarse- graining just a matter of 
taste or there arc physical or mathematical reasons to chose one or another? 
The reader must decide by himself. 

6.7 Appendix 6. A : Wigner function integral ([46], [4 7]). 

We have continuously jumped from the classical to the quantum case and 
back. Therefore it is interesting to present a theory to formalize these jumps 
and to make some applications of it. 

Let p be a density matrix of Liouvillc space JC — 7i x Ti and let {\q >} be 
the configuration or position basis of the Hilbert space H. The corresponding 
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Wigner function reads: 



pw{q,p) = TT-^ < g + \\p\q - A > e^'^^dX (6.^.1) 

It can be proved that: 

Lpw{q,p) =^~^ J <q + MLp\q - A > e^'^^dX + 0{h) (6.A2) 

where L is respectively the classical and quantum Liouville operator. In the 
classical limit ^ — > therefore pw can be considered as the classical distri- 
bution function corresponding to p. As in the classical regime we practically 
works in this limit we will consider that eq.(6.A.l) is the relation between 
the quantum density matrix and the classical distribution function. In fact, 
even if pw is not generally positive definite, using the Wigner integral from 
classical equation we can pas to quantum equation and vice versa, as a few 
examples will show. E.g., let us observe that: 

\\ Pw \h 1 1 Pw(q,p)dqdp^ ffi 4 

^JdqJ<q + X\p\q - A > 6{X)dX = trp ^ ^ 

so to the classical norm corresponds the quantum trace. Also: 

ipw\Ow) = 1 1 Pwiq,p)Ow{q,p)dqdp = 
T^~'^IdqJJ J < q + X\p\q- X >< q + p\0\q - p > . . . 

xe^'P^^+^'^dpdXdp = ..[^.^-V 
TT-^JdqJ < q + X\p\q - X >< q - X\0\q - X > dX^tr{pO) 

Therefore to the inner product in classical Liouville space corresponds the 
inner product in the quantum Liouville space. This fact complete the analogy 
between classical and quantum spaces implemented by the Wigner integral. 

As an exercise we can compute the classical distribution function cor- 
responding to density matrices pi{t) and P2{t) of eqs. (5.3.26) and (5.3.27). 
As these equations will be used in section 8 where we will use eq (6.5.2) to 
compute the entropy neglecting 0(A) we will do so in this exercise: 

Pm(t)=7r-i/<g + A|/o°°(pi^e-i(-i-)*|l><a;| . . 

+h.c.)duj\q - A > e^'P^dX ^ ' 
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where: 



< q\u >— 




(6.^1.6) 



where 2m — 1, Ti — 1, etc. After an easy calculation we obtain: 

Pwi{q,P, t) ^ Pi,(2p-y^)^e^^^(^-^')*e2^(^+^)« + h.c (6.^.7) 

We can see that the main values of this distribution function are obtained 
when p = y/tJi since for other values there are rapid oscillations. Analogously: 

pw2{t) ^n-^Pu J <q + X\l >< l\q-X> e^'^^dX ^ 6{p - ^)....{6.A.S) 

therefore also in this case all the effect is concentrated around the energy 
cui that will correspond, in the application of section 8, to the characteristic 
energy of nuclear reactions. 

7 Entropy in Curved Space Time. 

We have mention cosmology twice: 

i. -The cosmological event horizon could be a way to explain a universal 
graininess of nature. 

ii. -Fine graining time asymmetry is explained using a system with no 
exterior, namely the universe. 

Furthermore there is a cosmological arrow to investigate, so we cannot 
avoid cosmology in a complete discussion of our subject. 

Many years ago Mach thought us that most of the basic physical facts can 
only be explained only if we consider the universe as a whole, e. g.: if we want 
to explain why a system is a inertial one or not we must consider the whole 
universe, the system will be inertial if it is in uniform translatory motion 
with respect to the matter of the whole universe. The arrows of time are not 
exceptions, since they have a global nature. In fact, lumps are solved by the 
coffee in all places in the same time direction, here and in the Andromeda 
nebula. We must explain why it is so, and we will find the explanation only if 
we define the arrow of time in global cosmological models. Thus let us begin 
studying the notion of entropy in curved space-time because cosmological 
model are presented in this kind of spaces, [2]. 
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7.1 Thermodynamics in special relativity. 

For phenomenological reasons we can assume that the laws of thermodynam- 
ics are vahd in the special relativity proper system of coordinates S**^. From 
the relativity principle we then know that this laws are also valid in every 
inertial system S in translatory uniform motion with respect to S'^, provided 
the quantities involved in these laws would be transformed in a convenient 
way. In other word we would like to obtain the "Lorentz transformation" 
that makes invariant the following laws 
i.-The first law:: 



AE^AQ- AW. (7.1.1) 

where E is the energy, Q the heath and W the work, 
ii.- The second law: 

> ^ (7.1.2) 

where S is the entropy and T the temperature and the equality holds only 
for reversible evolution. To do this let us suppose that: 

i. - The pressure is isotropic namely it is equal in all directions and 

ii. - Let us temporarily use for simplicity axes chosen in such a way that 
the velocity u of the system S with respect to the system is parallel to 
the X-axis. 

Then from ordinary special relativity we know the coordinate transfor- 
mation equations for the following mechanical quantities.: 

i. -For the volume v: 

V = voVl - (7.1.3) 

ii. -For the pressure p: 

P = Po (7.1.4) 

iii. -For the energy E: 

E=^^ (7.1.5) 

Vl — 
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iv.-For the work W : 



dW = VT^dWo + d{Eo+poVo) (7.1.6) 

V 1 — 

where the quantities with subscript "0" refers to system S^. 

Then for the covariance of the first law eq. (7.1.1) it is necessary and 
sufficient that: 

Q = Vl^Qo (7.1.7) 

so we have obtained the transformation law of the heath. 

Let us now consider a thermic system at S^. We can accelerate this 
thermic system up to the velocity m in a reversible and adiabatic way, so the 
entropy of the system is not modified and we obtain: 

S = ^0 (7.1.8) 

Finally, from eqs. (7.1.7,8) it is evident that the second law, eq. (7.1.2) 
will be covariant iif: 

T = Vl - u'^To (7.1.9) 

so we have obtained the change of coordinate equation of all the basic ther- 
modynamical quantities. 

Let us now find the corresponding equations in four-dimensional language. 

The first law is just a form of the conservation of energy, therefore its 
four-dimensional version will be: 

Q^rpHU ^ Q (7.1.10) 

where T^^ is a convenient energy momentum tensor-(//, v, ... — 0, 1, 2, 3). 

To deduce the four-dimensional form of the second law let us consider 
a small volume of a thcrmodynamical fluid v and let us call the entropy 
density at the point where this element of volume is located in such a way 
that (f)v is the entropy of the element. If St is a small period of time the 
second law reads: 

d 5Q 



j^{ct>v)5t > -f (7.1.11) 



or (i,j... = 1,2,3): 
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where — Combining terms we have: 
but: 

^ = (7.1.14) 



and from eqs. (7.1.3) and (7.1.8): 



^0 



u 



2 



.(7.1.15) 



^ = ^ (7.1.16) 

T To ^ ^ 



Thus we obtain: 



where Sv — vSt is the coordinate four-dimensional volume element, which is 
equal to Svq — -^^=^y/T^^ii?6t — 5v (cf. eqs. (7.1.3), (7.1.14)), the proper 
four dimensional volume element, so we can use either one or the other. Thus 
if we define the flow of proper entropy or "entropy vector" as: 

we obtain the four-dimensional version of the second law: 

d.S'^Svo > ^ (7.1.19) 

-to 

which is valid for all inertial systems and we can put Sv instead of 5vo.- 

7.2 Thermodynamics in general relativity. 

Using the transcription rules to go from special relativity to general relativity, 
namely: 

Vnu 9nv, ^/x V^, 5v \f^5v, (7.2.1) 
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the first law, eq. (7.1.3) reads: 



V^T^'' = (7.2.2) 

or introducing the tensor density T'^" — ^/^T^'' we have the ordinary di- 
vergence: 

^^(T'^' + in = (7.2.3) 

where f^" is the pseudo tensor density of potential energy-momentum. This 
would be the general relativity covariant equation that shows the closest 
resemblance to eq. (7.1.3). 

Using the transcription rules on the second law, eq. (7.1.19) we obtain: 

W,S^V^5v > ^ (7.2.4) 

-to 

where, being all the factors scalars, we have, in fact, obtained a equation 
which is valid for all coordinate systems. Introducing the density = 
y/^S^, since V ^^,8^ — ■:^dny/^S^ — -^^d^S'^, this last equation reads: 

d.S'^Sv > ^ (7.2.5) 

-to 

which, again, is the general relativity covariant equation that shows the closet 
resemblance with the special relativity second law (7.1.19). 

Of course these are not the unique covariant generalization of the ther- 
modynamical laws of general relativity but they are the simplest and they 
lead to successful applications. 

7.3 Thermodynamics in cosmology. 

Let us consider a Robertson- Walker metric: 

ds^ = di^ + aW (7.3.1) 

where da is the comoving arc length and a the scale factor or the radius of 
the universe. 

If the energy momentum tensor corresponds to a isotropic fluid with den- 
sity Poo and pressure po the first law reads: 

j^{pooaHa) +poj^{,aHa) = (7.3.2) 
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where 6a is a comoving-coordinate three-dimensional volume. 

If we consider a comoving thermic fluid there will be not exchange of heath 
among the comoving volumes, and = ^ = (1, 0, 0, 0), so the second law, 
as expressed by eq. (7.2.5) reads: 

d^{<Pou''V^) > = ^(0oa') > 0....(7.3.3) 

where 0o is the proper entropy density and a the scale factor or radius of the 
universe. If we multiplies this equation by the constant coordinate comoving 
volume 5a we obtain: 

^((poa^Sa) > (7.3.4) 

at 

This equation gives the recipe to compute the entropy in that comoving frame 
of a Robertson- Walker metric: multiplied the local proper entropy density 
by the proper volume., which is, of course a very reasonable and natural 
result, that perhaps it is so natural it that can simply be assumed from the 
beginning, but now it is rigorously proved. 

Let us check this result with just one calculation: We now that in a 
radiation dominated universe temperature follows the law: 

T^To — (7.3.5) 

a 

that can be obtained integrating eq. (7.3.2) if we takepo = |poo ~ T'^, namely 
the radiation state equation, and that entropy of a black-body radiation in 
given by the formula: 

S = ^CsT^V. (7.3.6) 

where Cs is the Stefan coefficient, T the temperature and V the volume. If 
we substitute this two last equation into eq. (7.3.4) we see that the evolution 
of a radiation dominated universe is reversible, as can be expected. 

From these consideration we can deduce that the only effect produced 
by the expansion of the universe, in isotropic models, is the temperature 
decreasing. This is the only effect we must take into account below. 
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8 The Cosmological Problem. 



8.1 The problem of the time asymmetry. 

The problem of the existence of the arrows of time or, what is the same think, 
the time asymmetry of the universe, can be stated, as we explained in the 
introduction, in the following question: 

i.-How can it be that there is time asymmetry in the universe if all the 
relevant physical laws are time symmetric? 

In fact, universe has several time asymmetries, namely the various arrows 
of time: the thermodynamical one, the electromagnetical one, the psicolog- 
ical one, etc., while its main laws are time-symmetric (because, as usual in 
this kind of discussions, we will neglect the time-asymmetric laws of weak 
interactions [1], since it is very difficult to imagine a reason that explains 
macroscopic time asymmetry based only in the asymmetry of these laws) 

A second question is to explain the fact that all the arrows of time point 
in the same direction. 

In this section we would like to answer these questions, giving an adequate 
mathematical framework to the problem and using several, old and new, well 
known ideas ([2], [3], [4]) 

Let us first review the main equation of section 2. If the state of a physical 
system is described by p (being p classically the distribution function or 
quantum mechanically the density matrix) wc will call p^^'" =/Cp the state 
with reversed initial conditions (e.g.: if K is the Wigner operator of quantum 
mechanics then: ICp = i^pi^r'''.,[17],[15],[48]). We will say that the conditions 
at t — are time symmetric if p^''^'"{0) = /Cp(0) = p(0) and time asymmetric 
otherwise. If p{t) is the state of the universe at time t , the universe would 
have a time symmetric evolution with respect to i = if (cf. eq. (2.2.21)): 



But the universe has, in fact, a time- asymmetric evolution, at least with 
respect with some instant of time, that we call t — 0, such that: 



If the evolution equations, embodied in the universe liouvillian operator L, 
are time-symmetric, namely (cf. eq. (2.2.11)): 



]Cp{t) = p{-t) 



(8.1.1) 



)Cp{t) + p{-t) 




1CLK) = L 



(8.1.3) 
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to time symmetric conditions at t = will corresponds a time symmetric 
evolution (8.1.1) and to time asymmetric conditions will correspond time 
asymmetric evolutions like (8.1.2). In fact,: 



Pit) = e-^V(O) 



(8.1.4) 



therefore if the t — condition is time-symmetric we have: 



ICp{t) = €■ 



iKLKh 



/Cp(0) = e^";Cp(0) = p{-t) 



(8.1.5) 



since K, is an antilinear operator (namely ICi = —i). In the same way the time 
asymmetric case can be demonstrated. Then the observed time asymmetry 
of the universe evolution which obey eq, (8.1.2) can be explained only in two 
alternative ways: 

i.-Really cq (8.1.3) is not exact and there is a small, but relevant, time- 
asymmetric term in the liouvillian (e.g. cause perhaps by the weak interac- 



namely the initial state of the universe is not time symmetric. 

So, if we reject weak interactions, or any clever manipulation of the, 
otherwise time-symmetric physical laws, as the origin of time-asymmetry, 
we must necessarily consider eq. (8.1.6) as the only possible cause of this 
phenomenon. As,-in principle, asymmetry is a more generic property than 
symmetry (as complex numbers are more frequent than real ones) eq, (8.1.6) 
seems very natural and, therefore, this will be the idea that we will adopt in 
this section. If eq. (8.1.6) is valid, from eq. (8.1.5) we have: 



i.e. eq. (8.1.2), the equation we must prove. 

Finally, let us remark that the same explanation can be used to explain 
the other two fundamental asymmetries of nature V and C.In fact, if: 



tions) or: 



n.- 



/Cp(0) ^ p(0) 



(8.1.6) 



lCp{t) ^ e^^V(O) = p{-t) 




Pp(0)^p(0), 





we will have: 



Vp{t)^p{t), 
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even if: 

VLV^ = L, CLC^ = L (8.1.10) 

e.g.: Eq. (1.9,2) can be demonstrated if we postulate the existence of a small 
fluctuation between the amounts of matter and antimatter at the beginning 
of the universe. 

8.2 Entropy, Fluctuations, and Irreversibility. 

Let us first study the thermodynamical arrow of time. So, let us consider the 
entropy S as the state function that represent more eloquently the thermo- 
dynamical state of the universe {S can be computed using coarse-graining 
entropy, or fine-graining entropy). Wc know that the vast proportion of pos- 
sible states of the universe will be near the equilibrium state and will have 
the equilibrium entropy S*. Nevertheless we know that fluctuations around 
the equilibrium state, namely less probable unstable states near the equilib- 
rium, will spontaneously appear and we also know that, in these fluctuations 
states, entropy will be smaller than iS*. 

Anyhow steady equilibrium state satisfy Liouville equation: 

Lp^ = (8.2.1) 

For simplicity let us consider that there is just one equilibrium state in the 
universe, as it is very likely since the universe looks chaotic and,therefore,it 
is at least ergodic, therefore from eq. (8.1.3) we have: 

L/Cp, = ICLlOlCp, = (8.2.2) 

Therefore: 

/Cp* = p* (8.2.3) 

Thus if we have p(0) = p* we will have a time-symmetric evolution and no 
thermodynamical arrow of time (in fact, the universe will always remain in 
state p*). But for an unstable non-equilibrium-fluctuation state p ,in general, 
we will have- that /Cp ^ p. Therefore it is enough to assume that the universe 
began (al t — 0) in one of these states and we will have a time-asymmetric 
evolution and a thermodynamical arrow of time, because the initial entropy 
is S* < S*^, and therefore there will be growing of entropy, both to the past and 
to the future of t = 0, since entropy will try to reach the equilibrium entropy 
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in both directions, (in the exceptional case that the initial non-equilibrium 
unstable state would be such that /Cp(0) = p(0) at t = it would be fCp ^ p 
at a different time, close io t — that can, as well, be taken as the origin 
of time in eq, (8.1.7)). Then it is enough to suppose that the universe began 
in a non-equilibrium unstable state to obtain the thermodynamical arrow of 
time and the second law of thermodynamics, if we conventionally-consider 
only times t > (and conventionally call to this period the '' future" of t = 0). 

This low entropy initial state of the universe could be consider as a fluc- 
tuation. In fact, irregular fluctuation of the equilibrium entropy are present 
in systems with a finite number of particles [3], but vanish if this number 
goes to infinity. Then, these fiuctuations cannot be consider if we work with 
a distribution p in Liouville space, as we have done in these lectures, because 
these distributions are probabilities computed assuming an infinite number 
of particles (or an infinite number of copies of the system) . Fluctuations can 
be introduced in several way, e. g.: 

i. - Using Boltzmann entropy as in [7]. 

ii. - Working in a rigged space where distribution, corresponding to a finite 
number of particles, namely p's, built using a finite number of Dirac's deltas, 
can be consider, etc. 

We will not discuss this subject further here. 

About this solution, to the problem of the initial low entropy state of 
the universe, it can be argue that this initial fiuctuation is very unlikely 
[8], since the universe is very big, perhaps even infinite. Nevertheless, we 
shall prove, in the next section, that this conjecture is unnecessary, since the 
initial instability is naturally produced by the universe expansion, so really 
no fluctuations are needed, that is why we do not discuss fluctuations in these 
lectures. 

For isolated subsystems within the universe time asymmetry can be ob- 
tained in a similar way. In fact, we use to imagine that these subsystem 

(Gibbs ink drop spreading in a glass of water, or the perfume spreading into 
the room, etc., etc.) began in an unstable initial state with low entropy (a 
concentrated ink drop, all the perfume inside a bottle, etc., etc.). But these 
initial states are always produced, not by unlikely fluctuations but by exter- 
nal agencies (the ink or the perfume factories), that use energy to produce 
these concentrations that they obtain from other subsystems in unstable ini- 
tial states (chemical-unstable coal or nuclear-unstable isotopes, etc. etc.) 
that, in turn, obtain their energy, via a chain of unstable states (like those 
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of the stars), from the universe unstable initial state. Therefore, we conclude 
that all time-asymmetric processes have a cosmological origin. The only dif- 
ference is that in the case of a subsystem we have a reason to consider only 
times t > to, being the time to = the time of creation of the initial unsta- 
ble state of the subsystem, since time t < to corresponds to a period before 
the creation of the unstable states by the external agency (the concentration 
period of the ink drop or the bottle of perfume), where the subsystem is not 
isolated. The subsequent diffusion of the ink drop, the perfume, etc. will 
produce the growth of the thermodynamical entropy. 

The quantum analog of these reasonings can be found in papers [15] [49], 
[50] and [39]. 

The creation of a low entropy state is, therefore produce either by a ini- 
tial fluctuation, in the case of the universe (but we will see that the initial 
fluctuation is not necessary in the next subsection), or by an external agency, 
in the case of subsystems within the universe. Thus, neglecting for a mo- 
ment, the fluctuations, we will call a "conspiracy" to the appearance of a low 
entropy state not produced by an external agency. Then we can conclude 
that conspiracies do not exists in nature.In fact, let us consider a system in 
a low energy unstable state produced by external agencies (e.g. a glass store 
and an elephant). Any process, within the system, will produced a growth of 
entropy (e.g. the elephant moving by the store and braking all the glasses). 
This is an irreversible process. In fact, its time reverse process (a film of the 
motion of the elephant played backward) is full of conspiracies and therefore 
do not exists in nature. Irreversibility, therefore, can also be explained in 
this way in our formalism. 

8.3 The problem of the coordination of the arrows of 
time. 

Now we must solve the second problem that can be stated in the following 
question: 

a. -Why all the arrows of time point in the same direction? 
Also we would like to show that the initial fluctuation is not strictly 
necessary. 

To solve these problems we will consider that the cosmological arrow of 
time, namely the growth of the radius or scale factor of the universe a, is 
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the master arrow of time, that defines the direction of all the others. First 
we will show that the thermodynamical arrow of time, namely the tendency 
to obtain a final equilibrium, points in the same direction than the master 
arrow. 

Let S^.. be the equilibrium entropy and S{t) the actual entropy of the 
matter and radiation within the universe at time t (therefore now we will 
work with an open system since we exclude the entropy of the gravitational 
field) . The entropy gap: 

AS^S^-S{t) (8.3.1) 

would be minus the conditional entropy —Hc{p\p*) .according to eq. (3.6.2), 
in full agreement with general relativity, if we take into account the change 
of the universe temperature, as explained in section 7, namely: 




(8.3.2) 



where p{t,x) and p^{x) are the corresponding local distribution functions, X 
the phase space x & X a point of this space. The distribution functions are 
normalized as: 

j pdx = 1, j Pi,dx = 1, (8.3.3) 

Now we can consider that as in eq, (5.3.24): 

p{t) = + (pi + p2e- 2*)e-i* = + pAe-2* (8.3.4) 

where the second term of the r.h.s. is some kind of correction around the 
equilibrium term, with a dumping factor with a characteristic time Ri7~^. We 
will only consider the universe evolution after decoupling time, the universe 
will be matter dominated and 7"^ = tjsiRi will be the characteristic time 
of nuclear reactions, that make the matter within the star evolve toward 
thermal equilibrium with the cosmic microwave background. Eq. (8.3.4) can 
be considered only as a phenomenological equation, that can be obtained if 
we use coarse-graining technics and we neglect the Zeno and Khalfin effects; 
but we know there is a rigorous way to eliminate these effects, if we use the 
rigged Hilbert space formahsm (fine-graining technics) as in eq. (5.3.24). 
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Pa = Pi + P2e 2* is normalized as: 

/ p/^dx — / pdx — / p^,dx — 0.... (8.3.5) 

Jx Jx Jx 

This normalization is also a consequence of eq. (5.3.14'). 

We will consider that |pa| <S P*, or i » 7"^, namely that the fluctuation 
is small compared with the equilibrium distribution function. Then the en- 
tropy gap AS*, expanding the logarithm and neglecting unimportant terms, 
as in eq. (6. 5. 2), reads: 

2 

AS ^ e-^* / ^dx > (8.3.6) 

Jx 

Thus, when 7 = the growing of entropy variation disappear. To compute 
the time derivative of let us use the model of eqs. (5.3.31) and. (5.4.17), 
described at the end of subsection 5.3, then the last equation reads: 

Tf 



AS ^ e"^* / -^e^p\dx (8.3.7) 



where we have explicitated the time variation in the first exponential func- 
tion and in T'(i) being the rest of the quantities time constant, since we can 
neglect the second time variable term of pA with respect to the first con- 
stant one (or we could keep both terms with an small modification of the 
formulae). p\ and are independent of the temperature because they are 
related with the nuclear reaction processes only. From paper [31] (or eqs. 
(6. A. 7) and (6. A. 8)) we can introduce a reasonable simplification and sup- 
pose that the only important values of the last integral are those around 
the characteristic energy of the nuclear processes,then: 

AS^Ce~^T^^e^ (8.3.7) 

where C is a time independent constant. The temperature evolution will be 
dominated by the radiation within the universe and,therefore,will follow eq. 
(7.3.5) so: 

AS = C'e-^Vie^ (8.3.8) 

where C is another time independent constant. Now we can compute the 
time derivative that reads: 

Ai(t)=(-7-|^ + ^)A5 (8.3.9) 

Za J-odQ 
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where - = H(t) ^ t^^ is the Hubble coefficient. Since we are in the matter 
dominated period we have: 

a = ao (-] ' (8.3.10) 



Jo, 
thus: 

Eq. (8.3.8) shows two antagonic. effects (fig 5). The universe gravitational 
field, embodied in the positive coefficient (and in the term t~^), is the external 
agency that mostly try to take the system away from equilibrium, while, 
on the other hand, the nuclear reaction, embodied in 7 try to convey the 
system toward equilibrium (but the gravitational term try to establish 
equilibrium). This to effects are equal at a critical times tcrSuch that: 

7to+(^)=^f^)' (8.3.12) 



tcrJ 3To \t_ 

Usually this equation will have two positive roots tcri < t^^^ (fig. 6) 

It is premature to give physical numerical values to the parameters of 
the model. In fact, this model is extremely simplified, since it is based in 
an homogeneous space geometry while the decaying processes are produced 
within the stars, so what we really need is an inhomogeneous geometry to 
properly describe the phenomenon. However, with reasonable numerical val- 
ues (essentially taking uji » To, 7"^ to) we can obtain the following 
conclusions; 

a.-The ffist root is in the region t « so the ffist term of the l.h.s. of 

3To \ 5 



the last equation can be neglected to obtain tcri = to ^ ■ (This quantity, 
with minus sign, gives the third negative root) At this time the entropy gap 
has a minimum. 

b.-The second root is in the region t » toso the second term of the l.h.s. 

can be neglected to obtain tcr2 = ( ^srX^ ) ^o- At this time the entropy gap 
has a maximum. 

Then wc can state the following conclusions: 

i.-lf t < tcri then the second term of the l.h.s. of eq. (8.3.12) dominates 

A 5'< 0, and there is a big value for the entropy gap that is rapidly thermal- 
ized.. 
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ii. -If: tcri<t < tcr2, thcii AS> , the r.h.s. of eq. (8.3.12) dominates 
and there will be a growing of the entropy gap, produced by the universe 
expansion, that is driving the universe away from equilibrium. There will 
be a growing of complexity in this period, them particles, atoms, molecules, 
galaxies stars, planets, and living beings appear.. 

iii. -On the contrary, if t > tcr2 then AS'< , the first term of the l.h.s. of 
eq. (8.3.12) dominates ,the entropy gap will diminish and the universe goes 
toward its final equilibrium state, driven by the nuclear reaction processes, 
in agreement with paper [4]. All the structures within the universe decay 
and disappears. 

Therefore: 

limA,5= (8.3.13) 

Numerical estimations show that tcri << << tcr2, in such a way that the 
first period can be, some how, neglected since probably this period take place 
before decoupling time. (Also tcr2 3> to as in paper [51]) 

iv. -Eq, (8.3.8) shows how the universe expansion creates, in a continuous 
way, the universe instability and complexity. This fact make the initial fluc- 
tuation hypothesis unnecessary. This instability is created toward the future, 
defined as the direction of the universe expansion. Eq. (8.3.8) shows also 
how the local nuclear reaction try to restore equilibrium, in the same time 
direction. The thermodynamical arrow of time is the local tendency to ther- 
modynamical equilibrium (and not the total entropy growth, since our system 
is not isolated because the entropy of the gravitational field was not consid- 
ered) . Therefore the thermodynamical arrow coincide with the cosmological 
arrow. 

v. - All this reasonings are also valid before recombination time, where we 
must use a much bigger 7, because in that period we must consider reaction 

with much smaller characteristic time, in fact, smaller than recombination 
time. Since the period t < t^ri probably lays in this period perhaps the 
universe reach also a thermodynamical equilibrium, and we can use the ar- 
guments of reference [52] to show that the electromagnetical arrow of time 
coincide with the cosmological one. Also the dumping factor e~'^ can be 
obtained if we consider a pole in the lower half-plane of the unphysical sheet 
of the energy complex plane; thus we must use the upper rim of the positive 
real axis cut and, therefore, retarded solutions [35]. 
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vi. - Finally ourselves are just subsystem with unstable initial state, pro- 
duced by external agencies, like the ink drop or the bottle of perfume, there- 
fore our thermodynamical arrow, that can be identify with our psicological 
arrow, points also as the cosmological arrow. So all the arrows of time points 
in the same direction. 

vii. - Therefore wc have given a mathematical formalism to the answers of 
the to main questions about the universe time asymmetry. We believe that 
the presented solution is quite satisfactory, only much more physical exam- 
ples must be studied with the fine-graining method and some mathematical 
refinement are missing (like those of paper [48]). When these examples would 
be studied and this refinements will be implemented wc will have a definitive 
and rigorous answer to these, long standing, fundamental questions. 

9 Conclusions. 

After all these explanation and discussions we believe that we can drawn the 
following conclusions: 

i. -There are not compelling local-physical motivation to choose one tech- 
nic or the other. Therefore it is not easy to see how to find a local cross- 
experiment to settle the matter. Probably this cross-experiment not even 
exists, so really both technics are physically equivalent. 

ii. -Coarse-graining is more "physical", since it works directly in usual 
Hilbert space. The price to pay is the introduction of an object, which is 
really alien to the theory, the projector. This projector is essentially arbi- 
trary, so coarse grainig will not have a deep physical meaning until a natural 
graininess will be not find. 

iii. -Fine-graining is more "mathematic", since it works in Rigged Hilbert 
Space. But after paying this price, we are not force to introduce any object 
alien to the theory. In this sense fine-grainig is more pure, and really it cannot 
be distinguish from no-graining. Therefore it seems that fine-graining is con- 
ceptually superior even if, from the operational point of view, coarse-graining 
could be more convenient. Anyhow fine-graining has also its ambiguities, e.g.: 
the choice of the test function space, even if it seems more probable that we 
would find a canonical choice of this space, in the future, than a canonical 
choice of the coarse-graining projector. 

iv. For conceptually difficult chapters of physic, like cosmology or quan- 
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turn measurement theory, it is advisable to use fine- graining, since it is con- 
ceptually superior to coarse-graining. Thus, perhaps a global cross-experiment 
that shows the convenience of use one technic or the other could be find using 
cosmological reasonings. 
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11 Figures. 

0.-q{t) and p{t) functions for time-symmetric solutions, with respect tot — 0. 

1. -The baker transformation. 

2. -The P{t) graphic, showing Zeno effect, the exponential behavior, and, 
Khalfin effect. 

3. -The r curve- 

4. -The P' curve. 

5. -AjS' showing the minimum and the maximum. 

6. -A5" showing the two roots. 
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